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ABSTRACT 


Geochemical  reaction  path  modeling  is  useful  for  rapidly  assessing  the  extent  of 
water-aqueous-gas  interactions  both  in  natural  systems  and  in  industrial  processes.  Modeling  of 
some  systems,  such  as  those  at  low  temperature  with  relatively  high  hydrologic  flow  rates,  or 
those  perturbed  by  the  subsurface  injection  of  industrial  waste  such  as  C02  or  H2S,  must  account 
for  the  relatively  slow  kinetics  of  mineral-gas-water  interactions.  We  have  therefore  compiled 
parameters  conforming  to  a  general  Arrhenius-type  rate  equation,  for  over  70  minerals,  including 
phases  from  all  the  major  classes  of  silicates,  most  carbonates,  and  many  other  non-silicates.  The 
compiled  dissolution  rate  constants  range  from  -0.21  log  moles  in'"  s'  for  halite,  to  -17.44  log 
moles  m'  s-1  for  kyanite,  for  conditions  far  from  equilibrium,  at  25  °C,  and  pH  near  neutral. 
These  data  have  been  added  to  a  computer  code  that  simulates  an  infinitely  well-stirred  batch 
reactor,  allowing  computation  of  mass  transfer  as  a  function  of  time.  Actual  equilibration  rates 
are  expected  to  be  much  slower  than  those  predicted  by  the  selected  computer  code,  primarily 
because  actual  geochemical  processes  commonly  involve  flow  through  porous  or  fractured 
media,  wherein  the  development  of  concentration  gradients  in  the  aqueous  phase  near  mineral 
surfaces,  which  results  in  decreased  absolute  chemical  affinity  and  slower  reaction  rates.  Further 
differences  between  observed  and  computed  reaction  rates  may  occur  because  of  variables 
beyond  the  scope  of  most  geochemical  simulators,  such  as  variation  in  grain  size,  aquifer 
heterogeneity,  preferred  fluid  flow  paths,  primary  and  secondary  mineral  coatings,  and  secondary 
minerals  that  may  lead  to  decreased  porosity  and  clogged  pore  throats. 


1.  INTRODUCTION 


Equilibrium-based  geochemical  modeling  is  a  method  for  predicting  the  identity  and 
extent  of  chemical  reactions  in  both  geological  and  industrial  processes.  The  method  usually 
entails  computer-aided  simultaneous  solution  of  a  system  of  non-linear  equations.  Given  a 
pressure  (P),  a  temperature  (T),  and  a  bulk  chemical  composition  (x),  a  geochemical  modeling 
program  or  simulator  computes  the  equilibrium  distribution  of  the  chemical  components  among 
gaseous,  aqueous,  or  mineral  phases.  The  variables  P,  T,  or  x  can  be  changed  incrementally  over 
the  course  of  a  simulation  with  equilibrium  recomputed  at  each  step,  thereby  defining  a  reaction 
path.  Some  geochemical  simulators,  such  as  GAMSPATH  (for  a  description  of  GAMSPATH, 
and  comparison  of  GAMSPATH  to  other  geochemical  simulators,  see  Perkins  et  ah,  1997), 
require  the  presence  of  an  aqueous  fluid  phase  and  are  well  suited  to  simulation  of  systems  at 
shallow  to  moderate  depth  in  the  Earth’s  crust,  or  wherever  there  is  open  space  that  is  saturated 
with  water,  such  as  pores  or  fractures.  Water  is  the  medium  that  allows  mass  transfer  by 
diffusion  or  fluid  flow,  and  redistribution  of  the  components  via  chemical  reaction  among  the 
gaseous  and  various  mineral  phases,  and  species  in  the  aqueous  phase.  Considered  herein  is  the 
quantification  of  the  rates  at  which  these  reactions  occur. 

Local  non-redox  chemical  equilibria  among  fluids,  gases,  and  alteration  minerals  is  a 
valid  assumption  in  many  volcanic-hosted  hydrothermal  systems  (Arnorsson,  1983;  Arnorsson  et 
ah,  1983;  Giggenbach,  1980;  Giggenbach,  1981;  Pang  and  Reed,  1998;  Reed,  1982;  Reed,  1997; 


1 


Reed,  1998;  Reed  and  Spycher,  1984),  and  in  many  sedimentary  systems  with  low  fluid  flow 
rates  at  temperatures  as  low  as  75-80  °C  (Bazin  et  al.,  1997a;  Bazin  et  al.,  1997b;  Palandri  and 
Reed,  2001).  In  other  systems,  equilibrium  with  an  alteration  mineral  assemblage  cannot  always 
be  assumed,  especially  in  systems  where  low  temperature  leads  to  slow  rates  of  chemical 
reaction,  or  if  there  exists  a  perturbation  away  from  equilibrium,  e.g.  relatively  fast  hydrologic 
flow  rates  from  one  lithologic  regime  into  another.  It  is  important  to  recognize  that  equilibrium 
between  aqueous  fluid  and  primary  minerals  is  commonly  lacking,  and  in  some  cases  cannot  be 
attained,  e.g.  anorthite  or  forsterite  at  T  and  P  of  H?0  liquid-vapor  saturation.  Further, 
equilibrium  between  aqueous  fluid  and  secondary  minerals  containing  redox  elements  such  as  Fe 
and  S  generally  requires  higher  temperatures  than  for  those  minerals  containing  only  non-redox 
elements  (other  than  oxygen). 

For  many  geologic  systems  it  is  desirable  to  know  the  length  of  time  required  for  a 
system  to  equilibrate  with  respect  to  alteration  minerals,  and  perhaps  more  importantly,  the  rate 
of  dissolution  of  primary  minerals.  An  illustrative  example  is  the  geologic  sequestration  of 
carbon  dioxide  (C02)  from  anthropogenic  sources,  such  as  fossil  fuel-fired  electrical  power 
plants,  is  injected  into  depleted  oil  or  gas  reservoirs,  coal  seams,  or  deep  saline  aquifers.  As  a 
result  of  C02  injection,  these  systems  are  perturbed  far  from  equilibrium,  and  will  require  some 
length  of  time  to  re-equilibrate  with  alteration  minerals,  due  to  the  relatively  slow  kinetics  of 
mineral-C02-water  interactions  at  aquifer  temperatures  and  pressures.  While  dissolution  of 
supercritical  C02  into  aquifer  fluids  is  expected  to  be  kinetically  rapid,  the  rate  at  which 
dissolved  C02  reacts  with  primary  minerals,  and  the  rates  at  which  secondary  minerals 
precipitate  are  much  slower  (Gunter  et  al.,  1997). 

To  model  systems  where  the  time  dependence  of  geochemical  reactions  is  of  critical 
importance,  there  is  a  need  to  quantify  rates  of  mineral  dissolution  and  precipitation,  and  to  a 
lesser  extent,  of  aqueous  speciation.  Therefore,  we  compiled  directly  or  determined  using 
regression  and  curve-fitting,  the  parameters  conforming  to  a  general,  Arrhenius-type  rate 
equation  for  over  70  minerals,  including  most  of  the  rock-forming  minerals.  For  a  few  minerals, 
we  have  also  compiled  more  specialized  rate  equations  and  parameters,  where  some  workers 
have  clearly  shown  that  those  equations  provide  a  better  fit  to  experimental  data,  e.g.  quartz 
(Dove,  1994;  Dove,  1999;  Dove  and  Nix,  1997)  and  barite  (Dove  and  Czank,  1995).  The  rate 
parameters  compiled  herein  are  for  surface  controlled  reactions  at  conditions  that  are  far  from 
equilibrium,  rather  than  diffusion  controlled  reactions.  This  compilation  does  not  contain  rate 
parameters  for  aqueous  speciation. 

2.  METHODS 

2.1  Rate  Equations 


Any  rate  equation  that  is  to  be  used  in  geochemical  computer  modeling  must  contain  only 
parameters  that  are  available  during  program  execution  as  variables  or  constants.  The  general 
rate  equation  and  rate  parameters  compiled  herein  are  compiled  specifically  for  the  program 
GAMSPATH  (Perkins  et  al.,  1997),  which  uses  a  general,  semi-empirical  rate  equation  to  which 
experimental  rate  data  for  many  minerals  can  be  reasonably  well  fit.  A  general  fonn  of  the 
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equation  for  a  single  reaction  mechanism  modified  from  Lasaga  and  co-workers  (1984;  1995; 
1998)  is  given  by 


d^  =  -SA-Ae%\f(a,)g(AGr) 
at 


dm 

dt 


and  for  multiple  mechanisms 
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The  rate,  dm/dt,  is  in  units  of  mol  s'  .  The  surface  area  is  represented  by  SA  (m  ),  A  is  the 

Arrhenius  pre-exponential  factor  (mol  m'2  s'1),  E  is  the  activation  energy  (J  mof1),  T  is  the 

temperature  (K),  and  R  is  the  gas  constant. 

The  dimensionless  terms  ffflij )  are  functions  of  the  activity  of  the  ith  aqueous  species 

participating  in  the  jth  reaction  mechanism.  The  function  for  a  single  mechanism  is  commonly 

expressed  as  the  product  of  species  activities 


/(«<)  =  \\a" 

<  (3a) 

where  the  exponents  n;  represent  the  reaction  order  with  respect  to  species  a,  (Lasaga,  1995; 
Lasaga,  1998;  Lasaga  et  al.,  1994;  Oelkers  et  ah,  1994).  Lor  example,  the  expression 


r  s  \  -0.5  0.5 

f(ai)  =  aH+ape  3+  (3b) 

is  for  one  of  the  reaction  mechanisms  for  pyrite  (McKibben  and  Barnes,  1986),  discussed  further, 
below.  Positive  and  negative  exponents  denote  catalysts  and  inhibitors  respectively,  except  in 
the  case  of  hydroxide-catalyzed  mechanisms  (as  above  in  eqn.  (3b)),  which  are  quantified  in 
terms  of  oh+  raised  to  a  negative  exponent  (recall  that  for  aqueous  fluids,  aH+aoH-  =  10'14)  to 
simplify  calculations.  It  is  important  to  recognize  that  eqn.  (3b)  becomes  problematic  when  a,  is 
zero  and  /?,  is  positive,  because  the  rate  becomes  zero,  which  cannot  occur  for  some  mechanisms, 
or  if  a,  is  zero  and  /?,  is  negative,  because  a  numerical  singularity  occurs.  These  problems  can  be 
avoided  by  defining/with  adsorption  isotherms  (Blum  and  Lasaga,  1991;  Brady,  1992;  Casey 
and  Sposito,  1992;  Lurrer  and  Stumm,  1986;  Gautier  et  al.,  1994;  Kline  and  Logler,  1981; 
Pokrovsky  and  Schott,  1999;  Stillings  and  Brantley,  1995),  which  have  the  form 
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This  expression  is  an  example  is  an  isothenn  of  the  Langmuir  type  for  quartz  (Dove,  1994;  Dove, 
1999;  Dove  and  Nix,  1997),  where  KNa+  is  the  adsorption  coefficient  for  Na+,  and  kad  is  a 
constant  factor  by  which  the  overall  rate  constant  in  adjusted.  The  sign  of  kad  is  positive  for 
catalysts  and  negative  for  inhibitors.  GAMSPATH  currently  does  not  implement  adsorption 
isotherms,  and  uses  the  form  as  in  eqns.  (3a)  and  (3b). 

The  dimensionless  term  g(AGr)  is  a  function  of  the  chemical  affinity  to  account  for 
slowing  of  reaction  rates  as  equilibrium  is  approached,  and  may  be  expressed  as  (Lasaga,  1995; 
Lasaga,  1998;  Lasaga  et  ah,  1994;  Oelkers  et  ah,  1994): 


g(AG,,)  =  (l-fF'r  =(1- 


Omega  (Q  =  Q/K)  is  the  mineral  saturation  index  where  Q  is  the  activity  product,  and  K 
is  the  equilibrium  constant.  The  parameters  p,  and  q,-  are  empirical  and  dimensionless,  although 
Pi  can  be  predicted  from  transition  state  theory  if  the  details  of  the  reaction  mechanism  are 
known  (Lasaga,  1995;  Lasaga,  1998);  /?,■  and  qt  have  been  quantified  for  only  a  few  minerals 
(discussed  below),  and  for  those  minerals  only  for  a  single  mechanism;  in  all  other  cases,  the 
default  values  for  /;,  and  <y,-  are  to  a  first  approximation,  unity.  Uncertainties  in  values  for  pt  and 
q\  are  currently  the  source  of  much  uncertainty  in  the  length  of  time  to  system  equilibration 
computed  in  models  of  water-rock  interaction.  It  has  been  observed  for  some  minerals  that  there 
is  an  apparent  significant  slowing  of  the  dissolution  rate  at  low  degrees  of  undersaturation 
(Lasaga,  1995;  Lasaga,  1998;  Lasaga  et  ah,  1994).  For  example  gibbsite  (Nagy  and  Lasaga, 
1992)  and  albite  (Burch  et  ah,  1993),  show  a  stepped  shape  in  a  plot  of  dissolution  rate  versus 
chemical  affinity. 

The  overall  equation  implemented  in  GAMSPATH  is  a  sum  of  multiple  mechanisms 
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but  recast  in  terms  of  rate  constants  at  298.15  K  (25  °C)  and  pH  =  0  to  facilitate  comparison  of 
various  rates  among  minerals  at  near-surface  T 
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In  general,  the  most  well-studied  mechanisms  are  those  in  pure  H20  (neutral  pH),  and  those 
catalyzed  by  H+  (acid)  and  OH'  (base).  For  many  minerals,  the  full  equation  includes  a  term  for 
each  of  these  three  mechanisms: 
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The  OH'-catalyzed  mechanism  in  eqn.  (7)  uses  the  activity  of  H+  raised  to  a  reaction  order  with 
negative  sign,  to  simplify  data  reduction.  Additional  terms  are  added  to  or  removed  from  eqn. 

(7)  to  account  for  other  mechanisms  or  the  absence  of  data.  Other  mechanisms  include  those 
catalyzed  by  HCO3',  (especially  the  carbonate  minerals,  which  are  quantified  in  terms  of  either 
HCO3'  or  Pco2),  Fc  (pyrite,  pyrrhotite),  and  02  (pyrite).  In  the  case  of  pyrite  and  pyrrhotite,  the 
OH'-catalyzed  mechanism  is  apparently  slowed  by  presence  of  Fe3+.  This  effect  is  accounted  for 
by  substituting  a  term  in  the  fonn  of  eqn.  (3b)  for  the  H+  activity  and  reaction  order  in  the  third 
tenn  in  eqn.  (7).  It  should  be  recognized  that  although  the  pyrite  dissolution  rates  were  measured 
at  low  pH,  the  reaction  order  with  respect  to  H+  is  negative,  i.e.  the  reaction  is  OH'-catalyzed. 

Precipitation  rate  data  do  not  exist  for  most  minerals,  because  in  mineral  precipitation 
experiments,  undesired  metastable  reaction  products  usually  precipitate  instead  of  the  desired 
mineral,  especially  far  from  equilibrium  at  high  degrees  of  super-saturation.  One  example  of  this 
behavior  is  the  precipitation  of  amorphous  silica  from  fluids  that  are  grossly  supersaturated  with 
respect  to  quartz.  Further,  quartz  precipitation  is  prohibitively  slow  at  low  degrees  of  super¬ 
saturation.  Therefore,  we  use  the  principle  of  microscopic  reversibility  (Lasaga,  1998): 


K 


equilibrium 


^ dissolution 
^ precipitation 


K 
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(8) 


where  the  precipitation  rate  constant  can  be  calculated  from  the  dissolution  rate  constant  and  the 
equilibrium  constant.  However,  it  must  be  recognized  that  this  relation  requires  that  both  the 
forward  and  backward  reactions  must  both  proceed  by  a  single  reversible  mechanism.  Related  is 
the  phenomenon  of  nucleation,  where  some  degree  of  super-saturation  is  required  for 
precipitation  to  proceed  (Lebron  and  Suarez,  1996;  Nordeng  and  Sibley,  1994;  Normand  et  ah, 
2002;  Reddy,  1986;  Schoonen  and  Barnes,  1991b;  Shiraki  and  Brantley,  1995;  Sibley  et  al., 

1987;  Steefel  and  Van  Cappellen,  1990);  in  some  cases  this  can  be  substantial,  as  for  siderite 
(Greenberg  and  Tomson,  1992;  Jensen  et  al.,  2002).  Nucleation  may  be  either  homogeneous, 
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where  nucleation  occurs  spontaneously  in  a  liquid  phase,  or  heterogeneous,  where  nucleation  of 
a  new  phase  occurs  on  the  surface  of  a  pre-existing  solid  phase  (Lasaga,  1998).  For  most 
minerals,  the  degree  of  super-saturation  required  for  homogeneous  or  heterogeneous  nucleation 
has  not  been  quantified. 

For  further  general  discussions  of  mineral  dissolution  and  precipitation,  the  reader  may 
wish  to  consult  Anbeek  (1993),  Bertrand  et  al.  (1994),  Brady  and  Walther  (1989),  Casey  and 
Sposito  (1992),  Fleer  and  Johnston  (1986),  Gallup  (1998),  Gerard  et  al.  (1998),  Helgeson  (1970), 
Jeschke  and  Dreybrodt  (2002a),  Kline  and  Fogler  (1981),  Lasaga  (1984),  Lasaga  (1995),  Lasaga 
et  al.  (1994),  Oelkers  et  al.  (1994),  Petrovich  (1981a),  Petrovich  (1981b),  Schoonen  et  al.  (1998), 
Steefel  and  Lasaga  (1992),  Steefel  and  Van  Cappellen  (1990),  Sverdrup  (1990),  Talman  and 
Nesbitt  (1988),  Tuncay  et  al.  (2002),  Wieland  et  al.  (1988),  Wood  and  Walther  (1983),  Xiao  and 
Lasaga  (1994),  and  Xiao  and  Lasaga  (1996). 


2.2  Data  Reduction 


In  this  section,  we  outline  a  general  strategy  for  the  determination  of  rate  parameters 
conforming  to  eqn.  (6)  from  experimental  data,  using  the  minerals  albite  and  magnesite  as 
examples.  Albite  is  used  as  an  example  because  it’s  rate  data  clearly  conform  to  eqn.  (6)  in  one 
of  its  simplified  forms,  eqn.  (7),  where  rates  are  controlled  by  H+  and  OH'  catalysis  and 
temperature.  The  magnesite  example  illustrates  some  problems  that  may  arise  if  the  ability  to 
account  for  adsorption  isotherms  is  absent:  saturation  of  a  mineral  surface  with  a  (H+)  catalyst, 
and  the  affects  of  Na+,  which  may  catalyze  or  inhibit  dissolution  depending  on  pH. 

The  variables  that  have  the  greatest  effect  on  reaction  rates  are  temperature  and  pH,  and 
fortunately,  these  variables  are  almost  always  reported  in  by  experimentalists.  Many  minerals 
show  a  U-shaped  relation  (V-shaped  on  log  plots)  between  dissolution  rates  and  pH,  roughly  in 
parallel  with  their  solubilities.  For  a  given  temperature,  rates  are  generally  slowest  at  and  near 
the  pH  where  the  mineral  surface  has  a  net  charge  of  zero,  usually  within  one  or  two  pH  units  of 
neutral.  At  low  pH,  rates  increase  with  decreasing  pH  catalyzed  by  H+,  and  at  high  pH,  rates 
increase  with  increasing  pH,  catalyzed  by  OH'. 

An  example  of  this  behavior  is  shown  in  Figure  1  for  the  mineral  albite;  the  data  points  at 
100  to  300  °C  are  from  Heilman  (1994a),  and  at  25  °C,  from  Chou  and  Wollast  (1985).  The 
curves  in  Figure  1  are  obtained  from  the  data  (open  symbols)  using  unweighted  non-linear  curve 
fitting  (gray  curves)  and  unweighted  piecewise  linear  regression  (black  curves).  Non-linear 
curve  fitting  is  accomplished  by  fitting  all  of  the  data  to  a  modified  version  of  eqn.  (7): 
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(9) 


The  chemical  affinity  term  has  been  omitted,  because  we  have  selected  data  that  is  obtained  from 
experiments  perfonned  under  conditions  that  are  far  from  equilibrium  where  Q  approaches  zero, 
and  the  chemical  affinity  tenn  approaches  unity.  The  data  have  been  nonnalized  to  1.0  m  ,  and 
the  surface  area  term  is  omitted  as  well. 
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Piecewise  linear  regression  is  accomplished  by  first  segregating  the  data  into  acid, 
neutral,  and  basic  regions,  and  fitting  the  groups  of  data  separately  to  each  of  the  terms  in  eqn. 
(9).  For  albite  the  acidic,  neutral,  and  basic  regions  were  selected  by  Heilman  (1994a)  to  be  in 
the  ranges  of  pH  1.3  to  4.0,  5.6  to  8.2,  and  8.6  to  10.3  respectively.  For  the  acidic  and  basic 
mechanisms,  the  equation 


dm 

rate  =  — 
dt 


^ acid  ^ 


can  be  transformed  into  the  linear  relation 


(10) 


log  (rate)  =  log  (AacU)  -  E"^2.3025 RT )  ~  n«'pH  0  0 

and  for  the  neutral  mechanism,  the  following  relation  is  used 

log  {rate)  =  log  (Aneutral)  —  2025RT)  (12) 

The  black  dashed  straight  lines  in  Figure  1  show  the  results  for  each  of  the  three 
mechanisms  at  100°C,  and  the  black  curves  represent  their  sum.  It  can  be  seen  in  Figure  1  that 
non-linear  curve  fitting  provides  a  slightly  better  fit  to  the  data,  but  for  most  minerals,  lack  of 
agreement  in  the  data  among  different  workers  makes  this  difference  negligible. 

The  results  for  piecewise  linear  regression  and  non-linear  regression  for  albite  are  shown 
in  Table  1  both  in  terms  of  the  Arrhenius  pre-exponential  factor  A,  and  in  terms  of  the  rate 
constant  at  25  °C  and  pH  =  0.  Our  results  for  piecewise  linear  regression  are  in  agreement  with 
those  of  Heilman  (1994a)  and  Chou  and  Wollast  (1985).  We  recommend  application  of  these 
parameters  only  to  modeling  of  systems  within  the  range  of  conditions  under  which  the 
parameters  were  obtained,  i.e.  from  25  to  300  °C,  and  pH25oc  from  1.3  to  10.3.  Albite  is  one  of 
the  few  minerals  for  which  p  and  q  have  been  experimentally  determined,  but  only  in  alkaline 
fluids  (Alekseyev  et  al.,  1997)  at  pH  =  9;  these  parameters  are  also  shown  in  Table  1. 

A  second  illustrative  example  is  for  magnesite  dissolution  at  25  °C  (Pokrovsky  and 
Schott,  1999).  In  this  case,  magnesite  dissolution  rates  are  catalyzed  by  H+,  and  at  low  pH  the 
mineral  surfaces  become  saturated  with  H+  so  that  further  decrease  in  pH  does  not  cause  further 
increase  in  dissolution  rate.  Further,  dissolution  rates  are  apparently  catalyzed  or  inhibited  by 
Na'  depending  on  pH.  Finally,  magnesite  dissolution  rates  are  independent  of  OH'.  These 
experiments  were  conducted  at  pH  between  0. 19  to  10.13  with  pH  controlled  using  HC1,  NaOH, 
NaHC03,  and  Na2C03,  and  ionic  strength  controlled  by  NaCl.  The  results  are  shown  graphically 
in  Figure  2. 
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Figure  1 .  Rates  of  albite  dissolution  from  25  to  300  °C  and  pH  from  1 .3  to  10.3.  Data  at 
100,  200,  and  300  °C  (open  diamonds,  triangles  and  squares,  respectively)  from 
Heilman,  1994,  and  at  25  °C  (open  circles)  from  Chou  and  Wollast,  1985.  Dashed 
curves:  results  from  non-linear  regression.  Straight  black  lines:  results  from  piecewise 
linear  regression  of  100  °C  data.  Black  curves:  sum  of  results  from  piecewise  linear 
regression. 


Table  1.  Albite  Dissolution  Rate  Parameters. 


Acid  Mechanism  Neutral  Mechanism  Base  Mechanism 


aA 

blog  k 

CE 

dn 

aA 

blog  k 

CE 

aA 

blog  k 

CE 

dn 

Linear 

3.23E+01 

-9.87 

65.0 

0.457 

1.52E+00 

-12.04 

69.8 

2.96E-05 

-16.98 

71.0 

-0.572 

Regression 

Non-Linear 

1.66E+01 

-10.16 

65.0 

0.317 

6.66E-02 

-12.56 

65.0 

1.1  IE-04 

-15.60 

66.5 

-0.471 

Regression 

P 

q 

P 

q 

P 

q 

eChemical 

— 

— 

— 

— 

0.76 

90.0 

Affinity 

a.  Arrhenius  pre-exponential  factor  A,  mole  m  s  for  use  with  equation  (5). 

b.  Log  rate  constant  k  computed  from  A  and  E,  25°C,  pH=0,  log  moles  m~2  s'1  for  use  with  equation  (6). 

c.  Arrhenius  activation  energy  E,  kJ  mole  \ 

d.  Reaction  order  n  with  respect  to  H+. 

e.  Chemical  affinity  parameters  p  and  q  default  to  unity  if  not  specified. 
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In  the  pH  range  of  0  to  ~3,  there  is  only  a  minimal  dependence  of  rates  on  pH,  the 
reaction  order  with  respect  to  activity  of  H+  is  small,  and  increased  H+  activity  leads  to  only  a 
small  increase  in  rates  (Figure  2,  open  squares).  This  behavior  results  from  the  saturation  of 
reactive  surface  sites  with  H+,  and  can  be  described  using  adsorption  isotherms  (Pokrovsky  and 
Schott,  1999).  In  the  pH  range  of  3  to  5,  faster  rates  (Fig.  2,  open  diamonds)  apparently  result 
from  catalysis  by  Na+  at  ionic  strength  0.1  to  0.5,  as  compared  with  slower  rates  in  experiments 
with  low  Na+  concentration  (Fig.  2,  open  triangles)  at  ionic  strength  <0.1;  this  behavior  can  be 
described  using  adsorption  isothenns  for  Na+.  Between  pH  values  of  5  and  8,  there  is  minimal 
dependence  on  the  activity  of  H+  or  Na+  (Fig.  2,  open  triangles)  over  the  entire  range  of  ionic 
strength  from  0.01  to  0.5.  At  pH  greater  than  8,  Na+  apparently  inhibits  reaction  rates  (Fig.  2, 
open  circles)  at  ionic  strength  from  0.004  to  0.20);  again,  this  behavior  can  be  described  using 
adsorption  isotherms  for  Na+;  the  fastest  rates  in  this  pH  range  are  slightly  slower  than  the  rates 
in  the  near-neutral  pH  range  because  alkaline  pH  values  require  a  finite  amount  of  Na+  from 
NaOH. 

The  data  used  in  calculation  of  magnesite  rate  parameters  are  for  low  ionic  strength  only 
(<0. 1)  in  the  acidic  pH  range  of  3  to  5,  and  in  the  circum-neutral  pH  range  of  5-8  all  data  are 
used;  these  data  are  denoted  with  open  triangles  in  Figure  2.  For  non-linear  regression,  these 
data  are  simultaneously  fit  to  a  form  of  eqn.  (9)  from  which  the  tenn  for  the  basic  mechanism  has 
been  omitted;  the  data  were  obtained  from  experiments  at  25  °C  only,  so  the  activation  energy  is 
omitted  as  well 


dm 

dt 


/  £-298.15 K  m  £-298.15Zx 

V  acid  ^H+  neutral  ) 


(13) 


For  linear  regression,  the  data  are  fit  to  a  modified  form  of  eqn.  (11)  for  the  acidic  mechanism 

log  (rate)  =  log  KfJdl5K  -  nR+  pH  (14) 

and  for  the  neutral  mechanism,  a  simple  average  is  calculated,  yielding  the  rate  constants  at 
25  °C  and  an  acidic  mechanism  (positive)  reaction  order.  The  gray  curves  in  Figure  2  show  the 
results  of  the  non-linear  regression,  the  black  dashed  straight  lines  in  Figure  2  show  separately 
the  results  for  both  the  acidic  and  neutral  mechanisms,  and  the  black  curve  represents  their  sum. 

Since  data  are  lacking  for  the  temperature  dependence  of  magnesite  dissolution  rates,  the 
activation  energies  for  calcite  (discussed  below)  are  used  as  a  first  approximation,  allowing  for 
the  calculation  of  the  pre-exponential  factors.  The  chemical  affinity  factors  p  and  q  are  taken 
directly  from  Pokrovsky  and  Schott  (1999),  and  were  derived  from  experiments  in  which  the 
initial  solutions  contained  varying  amounts  of  Mg2+  and  CO32'.  Data  for  the  rate  parameters 
derived  from  the  regression  are  summarized  in  Table  2. 

Neither  piecewise  linear  regression  nor  non-linear  curve  fitting  provide  a  good  fit  to  the 
data  under  acidic  and  basic  pH  conditions,  because  of  the  inability  of  the  selected  rate  equations 
to  account  for  surface  adsorption  of  various  rate  catalyzing  and  inhibiting  species.  However,  the 
non-linear  curve  fitting  method  allows  the  calculation  of  rate  parameters  that  are  accurate  to  a 
first  approximation  over  a  wide  range  of  circum-neutral  pH. 
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Figure  2.  Rates  of  magnesite  dissolution  at  25  °C  and  pH  from  0.19  to  10.13.  Data 
25  °C  from  Pokrovsky  and  Schott,  1999;  open  triangles,  data  used  in  piecewise  linear 
regression;  open  squares,  other  data.  Dashed  lines,  results  from  piecewise  linear 
regression  of  acidic  and  neutral  mechanisms.  Black  curves,  sum  of  results  from 
piecewise  linear  regression.  Gray  curve,  results  from  non-linear  curve  fitting. 


Table  2.  Dissolution  Rate  Parameters. 


Acid  Mechanism 

Neutral  Mechanism 

aA 

blog  k 

CE 

dn 

aA 

blog  k 

CE 

Linear 

1.40E-04 

-6.38 

14.5 

0.589 

6.05E-06 

-9.34 

23.5 

Regression 

Non-Linear 

Regression 

1.62E+32 

32.21 

14.5 

0.448 

4.27E-10 

-9.37 

23.5 

P 

q 

P 

q 

eChemical 

Affinity 

4.00 

1.00 

-1  r  „ 

— 

b.  Log  rate  constant  k  computed  from  A  and  E,  25°C,  pH=0,  log  mole  m"2  s'1  for  use  with 
equation  (6). 

c.  Arrhenius  activation  energy  E,  kJ  mole"1. 

d.  Reaction  order  n  with  respect  to  H+. 

e.  Chemical  affinity  parameters  p  and  q  default  to  unity  if  not  specified. 
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2.3  Limitations  and  Uncertainties 


Application  of  kinetics  to  geochemical  modeling  is  subject  to  a  number  of  limitations 
resulting  from  uncertainties  in  the  experimental  methods  and  materials  from  which  the  rate 
parameters  are  derived,  and  from  the  modeling  software  used  for  the  simulations.  Ideally, 
experimental  data  would  be  available  over  wide  ranges  of  temperature  and  pressure  for 
dissolution  of  perfectly  crystalline  minerals  under  conditions  far  from  equilibrium  i.e.  where  the 
chemical  affinity  is  large,  so  that  further  increases  in  chemical  affinity  lead  to  no  observed 
increase  in  reaction  rates.  Additional  data  would  be  available  to  quantify  slowing  of  reaction 
rates  as  the  chemical  affinity  decreases,  as  equilibrium  is  approached  with  increasing 
concentration  of  dissolved  components  that  are  contained  in  the  mineral  structure.  Further  data 
would  be  available  for  dissolved  components  not  contained  in  the  mineral  structure  that  might 
catalyze  or  inhibit  the  reactions.  Further,  data  would  also  be  available  to  account  for  varying 
mineral  compositions  in  each  solid  solutions  series.  Finally,  data  would  be  available  for  minerals 
that  are  not  perfectly  crystalline,  i.e.  minerals  with  crystal  defects,  cation  disordering,  or  igneous 
exsolution  features  —  for  which  the  required  number  of  experiments  would  be  nearly  infinite. 
Most  of  these  data  do  not  yet  exist,  and  for  those  minerals  for  which  data  exist,  data  are  available 
mostly  for  experiments  conducted  only  under  conditions  that  are  far  from  equilibrium,  and 
values  for  p  and  q  are  assumed  to  be  equal  to  unity,  which  is  likely  incorrect. 

Additional  limitations  include  various  factors  of  experimental  design  and  solid  phase 
preparation:  1)  stirring  rates  may  be  insufficient  to  ensure  that  there  are  not  concentration 
gradients  near  mineral  surfaces,  which  may  slow  dissolution  reactions(Alkattan  et  ah,  1998; 

Metz  and  Ganor,  2001),  especially  for  minerals  that  have  fast  overall  rates  such  as  calcite 
(Alkattan  et  ah,  1998);  2)  in  some  experimental  reactor  designs,  particle  abrasion  may  lead  to 
faster  apparent  rates;  3)  grinding  of  samples  may  produce  surface  fines  which  can  be  removed 
with  solvents,  but  also  may  produce  particles  which  may  chemically  re-bind  to  mineral  surfaces 
(Petrovich,  1981b),  and/or  produce  crystal  defects  at  or  near  mineral  surfaces  (Bloom,  1983; 
Furrer  et  ah,  1993;  Nagy,  1995;  Nagy  et  ah,  1991;  Petrovich,  1981b;  Wieland  and  Stumm,  1992), 
which  may  lead  to  faster  initial  apparent  rates;  and  4)  there  is  also  a  degree  of  uncertainty  as  to 
whether  the  reactive  surface  area  can  be  equated  with  the  geometric  or  BET  surface  areas 
(Aagaard  and  Helgeson,  1982;  Gautier  et  ah,  2001;  Helgeson  et  ah,  1984;  Lasaga,  1998).  At 
least  one  study  (Gautier  et  ah,  2001)  of  quartz  dissolution  has  shown  that  an  increase  in  BET 
surface  area  during  dissolution  consists  primarily  of  un-reactive  etch  pit  walls,  and  concludes 
that  geometric  surface  area  provides  for  more  accurate  calculation  of  dissolution  rates. 

The  rate  equation  we  have  selected,  eqn.  (7),  is  also  subject  to  further  limitations.  It  does 
not  account  for  temperature  dependence  of  the  pre-exponential  factor  resulting  from  variation  in 
the  velocity  in  solution  of  dissolved  species,  nor  does  it  does  not  account  for  apparent  pH 
dependence  of  the  activation  energy,  which  can  be  described  using  adsorption  isotherms  (Lasaga, 
1995).  The  selected  equation  applies  only  to  bulk  dissolution  rates,  and  does  not  account  for 
crystallographic  control  of  mineral  dissolution,  i.e.  variation  of  dissolution  rates  on  different 
crystallographic  surface  planes. 
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3.  RESULTS 


3.1  Tectosilicates 

3.1.1  Si02  Polymorphs 

3. 1.1.1  Quartz 


Rate  data  and  parameters  for  quartz  from  several  sources  were  evaluated  (Bennett,  1991; 
Bennett  et  al.,  1988;  Bird  et  al.,  1986;  Blake  and  Walter,  1999;  Brady,  1992;  Dove,  1994;  Dove, 
1999;  Dove  and  Crerar,  1990;  Dove  and  Elston,  1992;  Dove  and  Nix,  1997;  Gautier  et  al.,  2001; 
Knauss  and  Copenhaver,  1995;  Knauss  and  Wolery,  1988;  Rimstidt  and  Barnes,  1980; 
Schwartzentruber  et  al.,  1987;  Tester  et  al.,  1994;  Welch  and  Ullman,  1992).  Data  from  Tester  et 
al.  were  selected  for  the  neutral  pH  mechanism  because  the  authors  were  able  to  correlate 
consistently  their  experimental  results  with  the  results  of  many  workers,  including  Bennett 
(1991),  Berger  et  al.  (1994),  Blum  et  al.  (1990),  Brady  and  Walther  (1990),  Dove  and  Crerar 
(1990),  Kitahara  (1960),  Rimstidt  and  Bames  (1980),  Siebert  et  al.  (1963),  Van  Lier  et  al.  (1960), 
and  Weill  and  Fyfe  (1964)  at  temperatures  from  25  to  625  °C.  In  the  acidic  region,  dissolution 
rates  are  apparently  independent  of  pH  (Dove,  1994;  Knauss  and  Wolery,  1988;  Rimstidt  and 
Barnes,  1980),  and  can  be  equated  with  rates  in  the  neutral  region  at  identical  T  and  P.  The  data 
used  to  regress  rate  parameters  using  eqn.  (1 1)  for  neutral  pH  are  shown  in  Table  3. 

Table  3.  Quartz  dissolution  Rate  Data  in  Pure  H2Q. 


T,  °C 

T,  K 

alog  k 

blog  k 

T,  °C 

T,  K 

alog  k 

blog  k 

23 

296.15 

-13.41 

-14.24 

400 

673.15 

-4.45 

— 

23 

296.15 

-13.38 

-14.21 

400 

673.15 

-4.35 

— 

50 

323.15 

-12.22 

-13.05 

440 

713.15 

-4.16 

— 

50 

323.15 

-12.21 

-13.04 

440 

713.15 

-4.35 

— 

50 

323.15 

-12.19 

-13.02 

440 

713.15 

-4.51 

— 

70 

343.15 

-11.33 

-12.18 

480 

753.15 

-4.41 

— 

100 

373.15 

-10.19 

-11.03 

480 

753.15 

-4.19 

— 

100 

373.15 

-10.22 

-11.06 

70 

343.15 

-12.16 

-12.22 

100 

373.15 

-10.20 

-11.04 

80 

353.15 

-11.46 

-11.52 

100 

373.15 

-10.45 

-11.3 

90 

363.15 

-11.16 

-11.22 

125 

398.15 

-9.32 

-10.17 

248 

521.15 

-6.08 

— 

125 

398.15 

-9.20 

-10.05 

250 

523.15 

-5.98 

— 

125 

398.15 

-9.08 

-9.93 

269 

542.15 

-5.84 

— 

125 

398.15 

-9.17 

-10.02 

291 

564.15 

-5.55 

— 

125 

398.15 

-9.19 

-10.03 

315 

588.15 

-5.22 

— 

125 

398.15 

-9.55 

-10.4 

332 

605.15 

-5.02 

— 

125 

398.15 

-9.75 

-10.6 

400 

673.15 

-3.95 

— 

125 

398.15 

-9.58 

-10.43 

625 

898.15 

-2.97 

— 

125 

398.15 

-9.82 

-10.67 

105 

378.15 

-10.27 

-11.5 

125 

398.15 

-9.59 

-10.44 

145 

418.15 

-8.36 

-9.59 

150 

423.15 

-8.42 

-9.26 

170 

443.15 

-8.43 

-9.66 

150 

423.15 

-8.48 

-9.33 

265 

538.15 

-6.56 

-7.79 

Table  3.  Quartz  dissolution  Rate  Data  in  Pure  H2Q  -  continued. 
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T,  °C 

T,  K 

alog  k 

blog  k 

T,  °C 

T,  K 

alog  k 

blog  k 

175 

448.15 

-8.38 

-9.23 

305 

578.15 

-6.14 

-7.37 

175 

448.15 

-8.34 

-9.19 

80 

353.15 

-10.75 

-11.18 

175 

448.15 

-8.38 

-9.23 

80 

353.15 

-10.89 

-11.25 

175 

448.15 

-8.28 

-9.12 

80 

353.15 

-10.99 

-11.52 

175 

448.15 

-8.32 

-9.16 

200 

473.15 

-7.63 

-7.93 

200 

473.15 

-7.88 

-8.72 

200 

473.15 

-7.35 

-7.65 

200 

473.15 

-7.82 

-8.67 

200 

473.15 

-7.28 

-7.58 

200 

473.15 

-7.05 

- 

200 

473.15 

-7.38 

-7.68 

200 

473.15 

-7.21 

- 

200 

473.15 

-7.26 

-7.56 

200 

473.15 

-7.30 

— 

201 

474.15 

-7.43 

-7.73 

250 

523.15 

-6.67 

- 

250 

523.15 

-6.74 

-7.04 

184 

457.15 

-7.96 

- 

300 

573.15 

-5.97 

-6.27 

184 

457.15 

-8.10 

- 

300 

573.15 

-6.01 

-6.31 

184.5 

457.65 

-8.32 

- 

300 

573.15 

-6.15 

-6.45 

202 

475.15 

-7.31 

— 

40 

313.15 

-13.13 

-12.7 

202.5 

475.65 

-7.49 

- 

50 

323.15 

-11.73 

-12.3 

203 

476.15 

-7.29 

— 

60 

333.15 

-11.33 

-11.9 

215 

488.15 

-7.05 

— 

70 

343.15 

-10.83 

-11.4 

221 

494.15 

-6.89 

— 

70 

343.15 

-10.93 

-11.5 

221 

494.15 

-6.96 

- 

200 

473.15 

-6.84 

-7.92 

221 

494.15 

-7.12 

- 

300 

573.15 

-5.38 

-6.46 

253 

526.15 

-6.38 

- 

300 

573.15 

-5.63 

-6.71 

255 

528.15 

-6.50 

- 

a.  Log  mole  m'2  s'1,  calculated  by  Tester  et  al.  (1994)  using  geometric  surface  area 

b.  Log  mole  m'2  s'1,  calculated  by  Tester  et  al.  (1994)  using  BET  surface  area 


Another  series  of  studies  from  Dove  and  co-workers  (Dove,  1994;  Dove,  1999;  Dove  and 
Crerar,  1990;  Dove  and  Elston,  1992;  Dove  and  Nix,  1997)  demonstrate  an  important 
phenomenon,  that  sodium  ion  catalyzes  the  OH'-catalyzed  (base)  dissolution  reaction  mechanism 
(as  do  Li,  K,  Mg,  Ca,  and  Ba).  These  data  were  not  used  because  adsorption  isotherms  must  be 
used  to  implement  their  results;  therefore,  we  use  the  results  from  Knauss  and  Wolery  (1988) 
directly  for  rates  under  alkaline  conditions,  which  are  given  terms  of  reaction  order  with  respect 
to  OH'. 

Results  from  the  regression  calculations  for  the  neutral  mechanism  using  both  geometric 
and  BET  surface  area  performed  herein  and  those  performed  by  Tester  et  al.  (1994),  and  results 
from  Knauss  and  Wolery  (1988)  for  the  base  mechanism  are  summarized  in  Table  4.  The  source 
of  the  small  difference  between  our  results  and  those  of  Tester  et  al.  (1994)  is  uncertain.  The 
chemical  affinity  parameters  were  not  detennined,  and  the  chemical  affinity  parameters  p  and  q 
are  assumed  to  be  unity. 
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Table  4.  Quartz  Dissolution  Rate  Parameters. 


Linear  Regression  Parameters 


Neutral  Mechanism  Base  Mechanism 


aA 

blog  k 

CE 

aA 

blog  k 

CE 

dn 

equartz 

333 

-13.40 

90.9 

- 

- 

- 

- 

fquartz 

276 

-13.34 

90.1 

~ 

~ 

~ 

~ 

9quartz 

23.3 

-13.99 

87.6 

~ 

~ 

~ 

~ 

hquartz 

24 

-13.99 

87.7 

- 

- 

- 

- 

'quartz 

— 

— 

— 

10 

-16.29 

108366 

-0.5 

a.  Arrhenius  pre-exponential  factor  A,  mole  m'2  s'1  for  use  with  equation  (5). 

b.  log  rate  constant  k  computed  from  A  and  E,  25°C,  pH=0,  log  mole  m'2  s'1  for  use  with  equation  (6). 

c.  Arrhenius  activation  energy  E,  kJ  mole'1. 

d.  Reaction  order  n  with  respect  to  H+. 

e.  Calculated  using  geometric  surface  area. 

f.  Calculated  by  Tester  et  al.  (1994)using  geometric  surface  area. 

g.  Calculated  using  BET  surface  area. 

h.  Calculated  by  Tester  et  al.  (1994)using  BET  surface  area. 

i.  From  Knauss  and  Wolery  (1988),  A  adjusted  here  from  491  to  10  to  be  consistent  with  the  results  of  Tester  et  al.  (1994). 


3.1 .1 .2  Amorphous  Si02,  Cristobalite,  and  Si02  Polymorph  Precipitation 

Rate  data  and  calculated  rate  equation  parameters  for  amorphous  silica  dissolution  from 
Icenhower  and  Dove  (2000),  and  Rimstidt  and  Bames  (1980)  were  evaluated.  Data  from 
Icenhower  and  Dove  (2000)  were  selected  because  they  conducted  many  individual  experimental 
runs,  and  they  were  able  to  correlate  their  data  with  data  from  other  workers.  The  rate  parameters 
are  summarized  in  Table  6.  The  presence  of  small  amounts  of  NaCl  in  solution  significantly 
enhances  dissolution  rates  (2 lx  at  0.05m  Incenhower  and  Dove,  2000),  this  can  be  accounted  for 
by  using  adsorption  isotherms,  but  cannot  be  described  with  our  selected  rate  equations  (eqns.  5 
and  6) 

Rate  parameters  for  a-cristobalite  and  P-cristobalite  are  compiled  directly  from  Rimstidt 
and  Barnes  (1980).  The  rate  parameters  were  calculated  by  the  authors  under  the  assumption 
that  the  precipitation  rates  for  all  of  the  SiCf  polymorphs  are  equal  (Rimstidt  and  Barnes,  1980), 
and  eqn.  (8).  These  data  are  also  summarized  in  Table  6. 


Table  5.  Amorphous  Silica  Dissolution  Rate  Data  in  Pure  H20. 


T,  °C 

T,  K 

alog  k 

T,  °C 

T,  K 

alog  k 

40 

313.15 

-11.36 

80 

353.15 

-10.47 

40 

313.15 

-11.17 

80 

353.15 

-10.47 

60 

333.15 

-10.75 

250 

523.15 

-6.59 

60 

333.15 

-10.84 

225 

498.15 

-6.86 

60 

333.15 

-10.63 

225 

498.15 

-6.96 

60 

333.15 

-10.57 

200 

473.15 

-7.4 

60 

333.15 

-10.65 

200 

473.15 

-7.34 

60 

333.15 

-10.99 

175 

448.15 

-7.89 

60 

333.15 

-10.97 

175 

448.15 

-7.74 
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Table  5.  Amorphous  Silica  Dissolution  Rate  Data  in  Pure  H2O  -  continued. 


T,  °C 

T,  K 

alog  k 

T,  °C 

T,  K 

60 

333.15 

-10.92 

150 

423.15 

-8.41 

60 

333.15 

-10.98 

150 

423.15 

-8.53 

60 

333.15 

-10.93 

225 

498.15 

-6.86 

60 

333.15 

-10.97 

225 

498.15 

-7.05 

60 

333.15 

-10.93 

200 

473.15 

-7.24 

60 

333.15 

-11.05 

200 

473.15 

-7.27 

60 

333.15 

-11.22 

175 

448.15 

-7.97 

60 

333.15 

-11.2 

175 

448.15 

-8.1 

60 

333.15 

-11.15 

150 

423.15 

-8.24 

60 

333.15 

-11.18 

125 

398.15 

-8.75 

60 

333.15 

-11.09 

100 

373.15 

-9.36 

80 

353.15 

-10.28 

100 

373.15 

-9.54 

a.  Log  mole  m"2  s'1  for  use  with  equation  (6);  data  for  experiments  in  de-ionized  water  only,  from 
Icenhower  and  Dove,  2000. 


Table  6.  Amorphous  Silica  and  Cristobalite  Dissolution,  and  Silica  Polymorph 
Precipitation  Rate  Parameters. _ 


Dissolution 

Neutral  Mechanism 

aA 

blog  k 

CE 

dAmorphous  Silica 

6.65E+00 

-12.23 

74.5 

eAmorphous  Silica 

1.01E+01 

-12.31 

76.0 

fAmorphous  Silica 

1.85E-01 

-12.77 

68.7 

fa-cristobalite 

1.20E-01 

-12.31 

65.0 

f(3-cristobalite 

7.60E-02 

-12.14 

62.9 

Precipitation 

fAII  Si02  polymorphs 

2.02E-01 

-9.42 

49.8 

a.  Arrhenius  pre-exponential  factor  A,  mole  m'2  s'1  for  use  with  equation  (5). 

b.  log  rate  constant  k  at  25°C  and  pH  =  0,  log  mole  m'2  s’1  for  use  with  equation  (6). 

c.  Arrhenius  activation  energy  E,  kJ  mole'1. 

d.  Reported  values  for  k  and  E  from  Incenhower  and  Dove  (2000),  A  calculated  herein. 

e.  A  and  E  regressed  using  data  from  Incenhower  and  Dove  (2000),  Table  5;  k  calculated  from  A  and  E  herein. 

f.  Reported  values  for  k  and  E  from  Rimstidt  and  Barnes  (1980);  A  calculated  from  k. 
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3.1.2  Feldspars 


3. 1.2.1  Plagioclase  Feldspars 

Plagioclase  dissolution  rate  data  were  evaluated  from  Alekseyev  et  al.  (1997),  Amrhein 
and  Suarez  (1988),  Amrhein  and  Suarez  (1992),  Berg  and  Banwart  (2000),  Bevan  and  Savage 
(1989),  Blake  and  Walter  (1996),  Blake  and  Walter  (1999),  Blum  and  Stillings  (1995),  Brady 
and  Walther  (1989),  Brantley  and  Stillings  (1994),  Burch  et  al.  (1993),  Busenberg  and  Clemency 
(1976),  Carroll  and  Knauss  (2001),  Casey  et  al.  (1991),  Casey  and  Sposito  (1992),  Chen  and 
Brantley  (1997),  Chen  et  al.  (2000),  Chou  and  Wollast  (1984),  Chou  and  Wollast  (1985),  Fenter 
et  al.  (2003),  Fleer  and  Johnston  (1986), Gautier  et  al.  (1994),  Hajash  et  al.  (1998),  Helgeson  et 
al.  (1984)  Holdren  and  Berner  (1979),  Holdren  and  Speyer  (1987),  Heilman  (1994b),  Heilman 
(1994a),  Heilman  (1995),  Knauss  and  Wolery  (1986),  Knauss  and  Copenhaver  (1995),  Lagache 
(1965),  Lagache  (1976),  Liittge  et  al.  (1999),  Manley  and  Evans  (1986),  Murakami  et  al.  (1998), 
Murphy  and  Helgeson  (1987;  1989),  Oelkers  et  al.  (1994),  Oelkers  and  Schott  (1995),  Oelkers 
and  Gislason  (2001),  Oelkers  (2001b),  Oxburgh  et  al.  (1994),  Rafal'skiy  et  al.  (1990),  Rafal'skiy 
and  Prisyagina  (1991),  Stillings  et  al.  (1995),  Stillings  et  al.  (1996),  Swoboda-Colberg  and 
Drever  (1993),  Talman  and  Nesbitt  (1988),  Taylor  et  al.  (2000a),  Tsuzuki  and  Suzuki  (1980), 
Teng  et  al.  (2001),  Utsunomiya  et  al.  (1999),  Welch  and  Ullman  (1992),  Welch  and  Ullman 
(1993),  Welch  and  Ullman  (1996),  and  Welch  and  Ullman  (2000). 

Unlike  the  plagioclase  endmember  albite  discussed  above,  dissolution  rates  for  the 
remainder  of  the  plagioclase  solid  solution  series  are  not  as  well  constrained.  Most  of  the 
available  data  are  biased  toward  low  temperature  and  near-neutral  to  acidic  pH  (Blum  and 
Stillings,  1995).  One  important  but  poorly  quantified  observation  is  that  dissolution  rates 
increase  continuously  with  increasing  anorthite  content.  This  effect  is  more  pronounced  for 
more  anorthitic  compositions;  at  25°C  and  pH  =  5,  rates  increase  by  approximately  one  order  of 
magnitude  from  albite  (log  k  ~  -12  log  mol  m"  s'  )  to  labradorite  (log  k  ~  -1 1  log  mol  m'  s'  ), 
and  by  almost  three  orders  of  magnitude  between  labradorite  and  anorthite  (log  k  ~  -8  log  mol 
m'  s'  )  (Blum  and  Stillings,  1995).  This  effect  diminishes  with  decreasing  pH,  and  the 
difference  in  the  logarithm  of  rate  approaches  a  linear  relation  with  respect  to  composition  at  pH 
=  2.  A  second  observation  that  has  not  been  well  quantified  is  that  dissolution  activation  energy 
decreases  with  increasing  anorthite  content. 

Our  general  strategy  is  to  identify  data  from  various  workers  that  are  consistent  with  one 
another  as  primary  constraints,  and  to  interpolate  pre-exponential  factors  and  activation  energies 
that  are  consistent  with  the  data.  Wherever  possible  rate  data  are  selected  from  experiments  in 
flow-through  reactors,  where  dissolution  rates  were  measured  using  the  concentration  of 
dissolved  silica  in  the  reactor  output  fluid,  and  where  silica  concentrations  were  allowed  to  reach 
steady-state,  thus  indicating  steady-state  thickness  of  any  possible  leached  layer.  One  issue  not 
addressed  in  the  cited  works  is  the  extent  to  which  cation  disorder  in  the  minerals  may  affect  the 
given  rates.  Data  were  selected  for  the  plagioclase  solid  solution  series,  and  are  summarized  in 
tables  as  follows:  albite,  Table  7;  oligoclase,  Table  8;  andesine,  Table  9;  labradorite,  Table  10; 
bytownite,  Table  11;  anorthite,  Table  12.  Primary  constraints  are  the  data  and  results  for  albite 
(T  =  25  to  300  °C,  pH  =  1.3  to  10.3)  and  anorthite  (T  =  8  to  95  °C,  pH  =  2.43  to  10.2),  with 
interpolation  obtained  graphically  for  the  components  of  the  solid  solution.  The  derived 
parameters  are  summarized  in  Table  13.  Chemical  affinity  parameters p  and  q  for  albite  are 
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0.760  and  90.0  respectively  (Alekseyev  et  al.,  1997),  but  their  use  in  modeling  should  be  limited 
to  conditions  near  the  experimental  conditions  under  which  they  were  obtained,  300  °C  and  pH  = 
9. 


Table  7.  Albite  Dissolution  Rate  Data. 


Acid  Mechanism 

Neutral  Mechanism 

Base  Mechanism 

T,  °C 

pH 

alog  k 

T,  °C 

pH 

log  k 

T,  °C 

pH 

alog  k 

25 

1.1 

-10.25 

25 

5.1 

-11.82 

25 

9.95 

-10.97 

25 

2.05 

-10.55 

25 

5.35 

-11.75 

25 

9.95 

-11.05 

25 

2.45 

-10.82 

25 

5.47 

-11.65 

25 

9.95 

-11.13 

25 

2.9 

-11.25 

25 

5.63 

-11.83 

25 

10.6 

-10.95 

25 

3 

-11.1 

25 

5.63 

-11.92 

25 

11.2 

-10.55 

25 

3.5 

-11.4 

25 

5.73 

-11.92 

25 

11.55 

-10.6 

25 

4.1 

-11.47 

25 

7.73 

-11.83 

25 

12.3 

-10.38 

100 

2.0 

-9.05 

25 

7.95 

-11.67 

100 

10.3 

-8.92 

100 

2.0 

-9.07 

100 

5.7 

-9.47 

200 

9.3 

-7.15 

100 

4.0 

-9.00 

100 

5.7 

-9.62 

300 

8.6 

-6.25 

100 

4.0 

-9.50 

100 

7.7 

-9.35 

300 

8.6 

-6.01 

100 

4.0 

-9.38 

100 

8.2 

-9.77 

300 

9.2 

-5.97 

200 

2.0 

-6.72 

172 

5.7 

-8.09 

300 

9.2 

-5.46 

200 

2.0 

-6.63 

191 

5.6 

-7.61 

300 

10.0 

-5.29 

200 

2.0 

-6.80 

200 

5.6 

-7.79 

200 

4.0 

-7.54 

200 

7.5 

-7.79 

300 

1.3 

-5.13 

222 

5.6 

-7.27 

300 

2.2 

-5.17 

300 

5.7 

-6.35 

300 

2.2 

-5.57 

300 

5.7 

-6.44 

300 

3.4 

-6.06 

300 

5.7 

-6.30 

300 

4.0 

-6.36 

300 

6.8 

-6.13 

300 

7.4 

-6.30 

300 

7.7 

-5.92 

300 

7.7 

-6.03 

a.  Log  mole  m  z  s'1.  Data  at  25  °C  from  Chou  and  Wollast  (1985),  for  Amelia  albite,  composition  not  specified.  Data  from  100  to  300 


°C  from  Heilman  (1984),  for  Amelia  Albite,  Ab10oAn0Or0  (1984). 
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Table  8.  Oligoclase  Dissolution  Rate  Data. 


an 

Composition 

ab 

or 

H 

O 

o 

bpH 

clogk 

Source 

20 

~ 

- 

25 

2 

-10.85 

Casey  et  al.,  1991 

20 

- 

- 

25 

2 

-10.44 

Casey  et  al.,  1991 

21 

- 

- 

25 

2 

-10.69 

Casey  et  al.,  1991 

13 

- 

- 

25 

3 

-11.47 

Oxburgh  et  al.,  1994 

20 

~ 

~ 

25 

3 

-10.94 

Holdren  and  Speyer,  1987 

13 

~ 

~ 

25 

3.6 

-11.70 

Oxburgh  et  al.,  1994 

13 

- 

- 

25 

4 

-11.92 

Oxburgh  et  al.,  1994 

13 

- 

- 

25 

4.6 

-12.22 

Oxburgh  et  al.,  1994 

13 

- 

- 

25 

5.1 

-12.10 

Oxburgh  et  al.,  1994 

24 

70 

6 

25 

5.3 

-11.59 

Busenberg  and  Clemency,  1976 

22 

73 

5 

25 

5.4 

-11.63 

Stillings  et  al.,  1996 

13 

- 

- 

25 

7 

-11.96 

Oxburgh  et  al.,  1994 

a.  -  =  Not  given. 

b.  Acid  mechanism,  pH  2-4.6;  neutral  mechanism,  pH  5.1-7. 

c.  Log  mole  m'2  s'1. 
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Table  9.  Andesine  Dissolution  Rate  Data. 


Composition 
an  ab  or 

T,  °C 

bpH 

clog  k 

Source 

46 

— 

— 

25 

2 

-10.16 

Casey  et  al.,  1991 

46 

— 

— 

25 

2 

-10.01 

Casey  et  al.,  1991 

46 

— 

— 

25 

3 

-10.64 

Holdren  and  Speyer,  1987 

49 

— 

— 

22 

3.1 

-9.78 

Welch  &  Ullman.,  1993 

49 

— 

— 

22 

3.1 

-9.27 

Welch  &  Ullman.,  1993 

46 

~ 

~ 

25 

3.1 

-10.65 

Oxburgh  et  al.,  1994 

47 

51 

2 

25 

3.13 

-10.62 

Stillings  et  al.,  1996 

49 

— 

— 

22 

3.8 

-10.12 

Welch  &  Ullman.,  1993 

49 

— 

— 

22 

4.1 

-9.81 

Welch  &  Ullman.,  1993 

46 

~ 

~ 

25 

4.1 

-11.15 

Oxburgh  et  al.,  1994 

47 

51 

2 

25 

4.11 

-11.15 

Stillings  et  al.,  1996 

47 

51 

2 

25 

5.08 

-11.24 

Stillings  et  al.,  1996 

46 

~ 

~ 

25 

5.1 

-11.24 

Oxburgh  et  al.,  1994 

49 

— 

— 

22 

5.3 

-10.47 

Welch  &  Ullman.,  1993 

43 

50 

8 

25 

5.3 

-11.86 

Busenberg  and  Clemency,  1976 

49 

— 

— 

22 

5.6 

-10.51 

Welch  &  Ullman.,  1993 

43 

50 

7 

25 

5.61 

-11.52 

Stillings  et  al.,  1996 

49 

— 

— 

22 

5.7 

-10.78 

Welch  &  Ullman.,  1993 

49 

— 

— 

22 

5.8 

-10.74 

Welch  &  Ullman.,  1993 

35 

65 

0 

250 

6.83 

-8.48 

Rafal'skiyand  Prisyagina,  1991 

47 

51 

2 

25 

7.22 

-11.21 

Stillings  et  al.,  1996 

46 

~ 

~ 

25 

7.3 

-11.29 

Oxburgh  et  al  1994 

35 

65 

0 

150 

7.46 

-9.39 

Rafal'skiyand  Prisyagina,  1991 

49 

— 

— 

22 

7.7 

-10.7 

Welch  &  Ullman.,  1993 

49 

— 

— 

22 

9 

-10.31 

Welch  &  Ullman.,  1993 

a.  -  =  Not  given. 

b.  Acid  mechanism,  pH  2-4.11;  Neutral  mechanism,  pH  5. 1-7. 7;  base  mechanism,  pH  9. 

c.  Log  mole  m'2  s'1. 
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Table  10.  Labradorite  Dissolution  Rate  Data. 


an 

Composition 

ab 

or 

T,  °C 

bPH 

clog  k 

Source 

51 

— 

— 

230 

0.41 

-6.44 

Tsuzuki  and  Suzuki,  1980 

51 

— 

— 

245 

0.45 

-6.48 

Tsuzuki  and  Suzuki,  1980 

51 

— 

— 

230 

0.81 

-6.35 

Tsuzuki  and  Suzuki,  1980 

51 

— 

— 

245 

0.84 

-6.71 

Tsuzuki  and  Suzuki,  1980 

51 

— 

— 

230 

0.96 

-6.71 

Tsuzuki  and  Suzuki,  1980 

51 

— 

— 

245 

0.98 

-6.71 

Tsuzuki  and  Suzuki,  1980 

51 

— 

— 

230 

1.21 

-7.6 

Tsuzuki  and  Suzuki,  1980 

51 

— 

— 

245 

1.23 

-7.48 

Tsuzuki  and  Suzuki,  1980 

51 

— 

— 

230 

1.67 

-7.48 

Tsuzuki  and  Suzuki,  1980 

51 

— 

— 

245 

1.69 

-7.44 

Tsuzuki  and  Suzuki,  1980 

60 

- 

-- 

25 

2 

-9.85 

Busenberg  and  Clemency,  1976 

61 

-- 

-- 

25 

3.2 

-10.6 

Casey  et  al.,  1991 

61 

-- 

-- 

25 

3.2 

-10.5 

Taylor  et  al.,  2000 

61 

~ 

~ 

25 

3.2 

-10.5 

Taylor  et  al.,  2000 

61 

-- 

~ 

25 

3.2 

-10.6 

Taylor  et  al.,  2000 

61 

~ 

~ 

25 

3.2 

-10.5 

Taylor  et  al.,  2000 

61 

-- 

~ 

25 

3.2 

-10.6 

Taylor 

et  al.,  2000 

60 

— 

— 

31 

3.2 

-9.45 

Carrol 

and  Knauss,  2001 

60 

— 

— 

60 

3.2 

-9.02 

Carrol 

and  Knauss,  2001 

60 

— 

— 

61 

3.2 

-8.37 

Carrol 

and  Knauss,  2001 

60 

— 

— 

100 

3.2 

-8.34 

Carrol 

and  Knauss,  2001 

60 

— 

— 

30 

3.21 

-9.53 

Carrol 

and  Knauss,  2001 

60 

— 

— 

130 

3.22 

-7.67 

Carrol 

and  Knauss,  2001 

60 

— 

— 

100 

3.24 

-8.11 

Carrol 

and  Knauss,  2001 

60 

— 

— 

130 

3.24 

-7.66 

Carrol 

and  Knauss,  2001 

60 

— 

— 

100 

3.25 

-8.09 

Carrol 

and  Knauss,  2001 

60 

— 

— 

60 

3.29 

-8.83 

Carrol 

and  Knauss,  2001 

60 

— 

— 

100 

3.29 

-8.31 

Carrol 

and  Knauss,  2001 

60 

— 

— 

130 

3.3 

-7.72 

Carrol 

and  Knauss,  2001 

60 

— 

— 

100 

3.37 

-8.52 

Carrol 

and  Knauss,  2001 

60 

— 

— 

150 

3.5 

-7.18 

Carrol 

and  Knauss,  2001 

57 

43 

0 

200 

4.88 

-8.48 

Lagache,  1965 

53 

44 

3 

25 

5.3 

-11.9 

Taylor  et  al.,  2000 

70 

— 

— 

80 

6 

-9.01 

Blake  and  Walter,  1999 

a.  -  =  Not  given. 

b.  Acid  mechanism,  pH  0.84-6. 

c.  Log  mole  m'2  s'1. 
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Table  1 1 .  Bytownite  Dissolution  Rate  Data. 


an 

Composition 

ab 

or 

T,  °C 

pH 

blog  k 

Source 

76 

- 

— 

25 

2 

-9.74 

Oxburgh  et  al.,  1994 

76 

- 

- 

25 

3 

-9.7 

Oxburgh  et  al.,  1994 

76 

- 

- 

25 

3.1 

-10.31 

Oxburgh  et  al.,  1994 

77 

- 

~ 

35 

3.2 

-11.13 

Welch  and  Ullman,  2000 

76 

- 

~ 

25 

3.8 

-10.91 

Oxburgh  et  al.,  1994 

77 

~ 

~ 

22 

4 

-8.69 

Welch  and  Ullman,  1993 

76 

- 

- 

25 

4.4 

-10.91 

Oxburgh  et  al.,  1994 

76 

- 

- 

25 

4.75 

-8.67 

Casey  et  al.,  1991 

77 

— 

— 

22 

4.9 

-8.89 

Welch  and  Ullman,  1993 

76 

~ 

~ 

25 

5 

-10.97 

Oxburgh  et  al.,  1994 

76 

~ 

~ 

25 

5.3 

-9.91 

Holdren  and  Speyer,  1987 

77 

- 

— 

5 

5.65 

-11.86 

Welch  and  Ullman,  2000 

77 

— 

— 

20 

5.65 

-11.4 

Welch  and  Ullman,  2000 

76 

24 

0 

25 

5.65 

-10.76 

Stillings  et  al.,  1996 

77 

22 

- 

25 

5.65 

-11.97 

Busenberg  and  Clemency,  1976 

77 

— 

— 

22 

5.7 

-10.2 

Welch  and  Ullman,  1993 

77 

— 

— 

22 

7.2 

-10.12 

Welch  and  Ullman,  1993 

a.  -  =  Not  given. 

b.  Acid  mechanism,  pH  2-5.3;  Neutral  mechanism,  pH  5.65-7.2. 

c.  Log  mole  nY2  s'1. 
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Table  12.  Anorthite  Dissolution  Rate  Data. 


an 

Composition 

ab 

or 

T,  °C 

pH 

blog  k 

Source 

100 

0 

0 

25 

2 

-10.32 

Fleer,  1982 

100 

0 

0 

50 

2 

-9.72 

Fleer,  1982 

100 

0 

0 

70 

2 

-9.59 

Fleer,  1982 

92 

— 

— 

25 

3 

-9.77 

Holdren  and  Speyer,  1987 

93 

— 

— 

25 

3.6 

-10.24 

Amrhein  and  Suarez,  1992 

unk 

— 

— 

25 

4 

-10.65 

Bailey,  1974 

unk 

— 

— 

45 

4 

-10.35 

Bailey,  1974 

unk 

— 

— 

65 

4 

-10.35 

Bailey,  1974 

93 

— 

— 

25 

4.45 

-10.83 

Amrhein  and  Suarez,  1992 

93 

— 

— 

25 

4.65 

-11.33 

Amrhein  and  Suarez,  1992 

93 

— 

— 

25 

4.7 

-11.27 

Amrhein  and  Suarez,  1992 

93 

— 

— 

25 

6 

-10.14 

Amrhein  and  Suarez,  1992 

93 

— 

— 

25 

6 

-10.69 

Amrhein  and  Suarez,  1992 

100 

0 

0 

25 

6.7 

-10.83 

Berg  and  Banwart,  2000 

100 

0 

0 

25 

7.4 

-11.15 

Berg  and  Banwart,  2000 

93 

— 

— 

25 

7.65 

-11.37 

Amrhein  and  Suarez,  1992 

100 

0 

0 

25 

8.1 

-11.11 

Berg  and  Banwart,  2000 

100 

0 

0 

25 

8.4 

-10.82 

Berg  and  Banwart,  2000 

100 

0 

0 

70 

10.2 

-11.41 

Fleer,  1982 

97 

— 

— 

45 

2.46 

-6.79 

Oelkers  and  Schott,  1995 

97 

— 

— 

45 

2.47 

-6.86 

Oelkers  and  Schott,  1995 

97 

— 

— 

45 

2.65 

-7.09 

Oelkers  and  Schott,  1995 

97 

— 

— 

45 

2.56 

-7.01 

Oelkers  and  Schott,  1995 

97 

— 

— 

45 

2.56 

-7.06 

Oelkers  and  Schott,  1995 

97 

— 

— 

45 

2.58 

-6.97 

Oelkers  and  Schott,  1995 

97 

— 

— 

45 

2.48 

-6.94 

Oelkers  and  Schott,  1995 

97 

— 

— 

60 

2.53 

-6.65 

Oelkers  and  Schott,  1995 

97 

— 

— 

60 

2.78 

-7.09 

Oelkers  and  Schott,  1995 

97 

— 

— 

60 

2.54 

-6.75 

Oelkers  and  Schott,  1995 

97 

— 

— 

60 

2.65 

-6.94 

Oelkers  and  Schott,  1995 

97 

— 

— 

60 

2.56 

-6.82 

Oelkers  and  Schott,  1995 

97 

— 

— 

60 

2.52 

-6.76 

Oelkers  and  Schott,  1995 

97 

— 

— 

60 

2.52 

-6.74 

Oelkers  and  Schott,  1995 

97 

— 

— 

60 

2.51 

-6.58 

Oelkers  and  Schott,  1995 

97 

— 

— 

60 

2.72 

-6.92 

Oelkers  and  Schott,  1995 

97 

— 

— 

60 

2.54 

-6.77 

Oelkers  and  Schott,  1995 

97 

— 

— 

60 

2.77 

-7.04 

Oelkers  and  Schott,  1995 

97 

— 

— 

60 

2.52 

-6.70 

Oelkers  and  Schott,  1995 

97 

— 

— 

60 

2.62 

-6.87 

Oelkers  and  Schott,  1995 

97 

— 

— 

60 

2.50 

-6.68 

Oelkers  and  Schott,  1995 

97 

— 

— 

60 

2.65 

-6.96 

Oelkers  and  Schott,  1995 

97 

— 

— 

60 

2.78 

-7.16 

Oelkers  and  Schott,  1995 

97 

— 

— 

60 

2.53 

-6.78 

Oelkers  and  Schott,  1995 

22 


Table  12.  Anorthite  Dissolution  Rate  Data  -  continued. 


an 

Composition 

ab 

or 

H 

O 

o 

pH 

blog  k 

Source 

97 

— 

— 

60 

2.74 

-7.13 

Oelkers  and 

Schott,  1 995 

97 

— 

— 

60 

2.68 

-7.00 

Oelkers  and 

Schott,  1995 

97 

— 

— 

60 

2.63 

-6.94 

Oelkers  and 

Schott,  1 995 

97 

— 

— 

60 

2.81 

-7.21 

Oelkers  and 

Schott,  1995 

97 

— 

— 

60 

2.59 

-6.93 

Oelkers  and 

Schott,  1 995 

97 

— 

— 

60 

2.56 

-6.89 

Oelkers  and 

Schott,  1995 

97 

— 

— 

60 

2.71 

-7.14 

Oelkers  and 

Schott,  1 995 

97 

— 

— 

60 

2.54 

-6.90 

Oelkers  and 

Schott,  1995 

97 

— 

— 

60 

2.56 

-6.64 

Oelkers  and 

Schott,  1995 

97 

— 

— 

60 

2.60 

-6.77 

Oelkers  and 

Schott,  1995 

97 

— 

— 

60 

2.75 

-7.06 

Oelkers  and 

Schott,  1 995 

97 

— 

— 

60 

2.66 

-6.92 

Oelkers  and 

Schott,  1 995 

97 

— 

— 

60 

2.51 

-6.70 

Oelkers  and 

Schott,  1995 

97 

— 

— 

60 

2.55 

-6.85 

Oelkers  and 

Schott,  1 995 

97 

— 

— 

60 

2.70 

-7.13 

Oelkers  and 

Schott,  1 995 

97 

— 

— 

60 

2.61 

-7.02 

Oelkers  and 

Schott,  1995 

97 

— 

— 

60 

2.72 

-6.96 

Oelkers  and 

Schott,  1 995 

97 

— 

— 

60 

2.58 

-6.72 

Oelkers  and 

Schott,  1 995 

97 

— 

— 

60 

2.77 

-7.10 

Oelkers  and 

Schott,  1 995 

97 

— 

— 

60 

2.65 

-6.89 

Oelkers  and 

Schott,  1 995 

97 

— 

— 

60 

2.63 

-6.85 

Oelkers  and 

Schott,  1 995 

97 

— 

— 

60 

2.46 

-6.70 

Oelkers  and 

Schott,  1 995 

97 

— 

— 

60 

2.41 

-6.55 

Oelkers  and 

Schott,  1 995 

97 

— 

— 

60 

2.62 

-6.81 

Oelkers  and 

Schott,  1 995 

97 

— 

— 

60 

2.51 

-6.61 

Oelkers  and 

Schott,  1 995 

97 

— 

— 

60 

2.55 

-6.72 

Oelkers  and 

Schott,  1 995 

97 

— 

— 

60 

2.61 

-6.88 

Oelkers  and 

Schott,  1 995 

97 

— 

— 

60 

2.53 

-6.74 

Oelkers  and 

Schott,  1 995 

97 

— 

— 

60 

2.58 

-6.94 

Oelkers  and 

Schott,  1 995 

97 

— 

— 

60 

2.67 

-6.77 

Oelkers  and 

Schott,  1 995 

97 

— 

— 

60 

3.08 

-7.52 

Oelkers  and 

Schott,  1 995 

97 

— 

— 

60 

2.90 

-7.22 

Oelkers  and 

Schott,  1 995 

97 

— 

— 

60 

2.87 

-7.16 

Oelkers  and 

Schott,  1 995 

97 

— 

— 

60 

3.00 

-7.40 

Oelkers  and 

Schott,  1 995 

97 

— 

— 

60 

2.91 

-7.28 

Oelkers  and 

Schott,  1 995 

97 

— 

— 

60 

2.94 

-7.38 

Oelkers  and 

Schott,  1 995 

97 

— 

— 

75 

2.52 

-6.53 

Oelkers  and 

Schott,  1 995 

97 

— 

— 

75 

2.54 

-6.63 

Oelkers  and 

Schott,  1 995 

97 

— 

— 

75 

2.71 

-6.94 

Oelkers  and 

Schott,  1 995 

97 

— 

— 

75 

2.56 

-6.69 

Oelkers  and 

Schott,  1995 

97 

— 

— 

75 

2.61 

-6.78 

Oelkers  and 

Schott,  1995 

97 

— 

— 

75 

2.67 

-6.86 

Oelkers  and 

Schott,  1 995 
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Table  12.  Anorthite  Dissolution  Rate  Data  -  continued. 


an 

Composition 

ab  or 

T,  °C 

pH 

blog  k 

Source 

97 

— 

75 

2.51 

-6.63 

Oelkers  and  Schott,  1995 

97 

.. 

75 

2.50 

-6.61 

Oelkers  and  Schott,  1995 

97 

.. 

95 

2.65 

-6.71 

Oelkers  and  Schott,  1995 

97 

.. 

95 

2.60 

-6.61 

Oelkers  and  Schott,  1995 

97 

.. 

95 

2.53 

-6.48 

Oelkers  and  Schott,  1995 

97 

.. 

95 

2.73 

-6.91 

Oelkers  and  Schott,  1995 

97 

.. 

95 

2.58 

-6.59 

Oelkers  and  Schott,  1995 

97 

-- 

95 

2.55 

-6.54 

Oelkers  and  Schott,  1995 

97 

-- 

95 

2.68 

-6.80 

Oelkers  and  Schott,  1995 

97 

.. 

95 

2.58 

-6.64 

Oelkers  and  Schott,  1995 

Table  13.  Plagioclase  Dissolution  Rate  Parameters. 


Acid  Mechanism 

Neutral  Mechanism 

Base  Mechanism 

alog  k 

bE 

cn 

alog  k 

bE 

alog  k 

bE  cn 

Albite 

-10.16 

65.0 

0.457 

-12.56 

69.8 

-15.60 

71.0  -0.572 

Oligoclase 

-9.67 

65.0 

0.457 

-11.84 

69.8 

- 

- 

Andesine 

-8.88 

53.5 

0.541 

-11.47 

57.4 

- 

- 

Labradorite 

-7.87 

42.1 

0.626 

-10.91 

45.2 

- 

.. 

Bytownite 

-5.85 

29.3 

1.018 

-9.82 

31.5 

- 

-- 

Anorthite 

-3.50 

16.6 

1.411 

-9.12 

17.8 

- 

-- 

a.  Rate  constant  k  computed  from  A  and  E,  25°C,  pH  =  0,  mole  m'2  s 1 

b.  Arrhenius  activation  energy  E,  kJ  mole"1. 

c.  Reaction  order  n  with  respect  to  H+ 


3. 1.2.2  K-feldspar 

Dissolution  rate  data  were  evaluated  from  Alekseyev  et  al.  (1997),  Bevan  and  Savage 
(1989),  Blake  and  Walter  (1996),  Blake  and  Walter  (1999),  Blum  and  Stillings  (1995), 
Busenberg  and  Clemency  (1976),  Fenter  et  al.  (2003),  Gautier  et  al.  (1994),  Helgeson  et  al. 
(1984),  Holdren  and  Speyer  (1987),  Knauss  and  Copenhaver  (1995),  Lagache  (1976),  Manley 
and  Evans  (1986),  Rafal'skiy  et  al.  (1990),  and  Teng  et  al.  (2001).  Rate  parameters  are  taken 
directly  from  Blum  and  Stillings  (1995)  for  the  acidic  mechanism,  and  are  regressed  herein  for 
the  neutral  mechanism  from  data  (Table  14)  in  Helgeson  et  al.  (1984)  and  for  the  basic 
mechanism  from  data  (Table  14)  in  Bevan  and  Savage  (1989),  Blum  and  Stillings  (1995), 
Gautier  et  al.  (1994),  Knauss  &  Copenhaver  (1995).  The  derived  rate  parameters  are 
summarized  in  Table  15. 
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Table  14.  K-feldspar  Dissolution  rate  Data. 


an 

Composition 

ab 

or 

T,  °C 

bPH 

clog  k 

Source 

0 

16 

84 

25 

4 

-8.65 

Helgeson  et  a 

.  1984 

0 

16 

84 

50 

4 

-7.54 

Helgeson  et  a 

.  1984 

0 

16 

84 

100 

4 

-5.76 

Helgeson  et  a 

.  1984 

0 

16 

84 

150 

4 

-4.39 

Helgeson  et  a 

.  1984 

0 

16 

84 

200 

4 

-3.3 

Helgeson  et  a 

.  1984 

0 

16 

84 

250 

4 

-2.42 

Helgeson  et  a 

.  1984 

0 

16 

84 

300 

4 

-1.69 

Helgeson  et  a 

.  1984 

0 

16 

84 

350 

4 

-1.07 

Helgeson  et  a 

.  1984 

0 

16 

84 

100 

4 

-6.2 

Helgeson  et  a 

.  1984 

0 

16 

84 

150 

4 

-4.38 

Helgeson  et  a 

.  1984 

0 

16 

84 

200 

4 

-3.52 

Helgeson  et  a 

.  1984 

0 

16 

84 

200 

4 

-3.27 

Helgeson  et  a 

.  1984 

0 

16 

84 

200 

4 

-3 

Helgeson  et  a 

.  1984 

0 

16 

84 

250 

5.58 

-8.61 

Helgeson  et  a 

.  1984 

0 

16 

84 

250 

5.58 

-2.42 

Helgeson  et  a 

.  1984 

0 

16 

84 

200 

5.64 

-9.03 

Helgeson  et  a 

.  1984 

0 

16 

84 

200 

5.64 

-3.3 

Helgeson  et  a 

.  1984 

0 

16 

84 

300 

5.65 

-8.26 

Helgeson  et  a 

.  1984 

0 

16 

84 

300 

5.65 

-1.69 

Helgeson  et  a 

.  1984 

0 

16 

84 

150 

5.82 

-9.54 

Helgeson  et  a 

.  1984 

0 

16 

84 

150 

5.82 

-4.39 

Helgeson  et  a 

.  1984 

0 

16 

84 

350 

5.91 

-7.97 

Helgeson  et  a 

.  1984 

0 

16 

84 

350 

5.91 

-1.07 

Helgeson  et  a 

.  1984 

0 

16 

84 

25 

6 

-8.65 

Helgeson  et  a 

.  1984 

0 

16 

84 

50 

6 

-7.54 

Helgeson  et  a 

.  1984 

0 

16 

84 

100 

6 

-5.76 

Helgeson  et  a 

.  1984 

0 

16 

84 

150 

6 

-4.39 

Helgeson  et  a 

.  1984 

0 

16 

84 

200 

6 

-3.3 

Helgeson  et  a 

.  1984 

0 

16 

84 

250 

6 

-2.42 

Helgeson  et  a 

.  1984 

0 

16 

84 

300 

6 

-1.69 

Helgeson  et  a 

.  1984 

0 

16 

84 

350 

6 

-1.07 

Helgeson  et  a 

.  1984 

0 

16 

84 

100 

6.13 

-10.18 

Helgeson  et  a 

.  1984 

0 

16 

84 

100 

6.13 

-5.76 

Helgeson  et  a 

.  1984 

0 

16 

84 

50 

6.64 

-11.01 

Helgeson  et  a 

.  1984 

0 

16 

84 

50 

6.64 

-7.54 

Helgeson  et  a 

.  1984 

0 

16 

84 

25 

7 

-11.52 

Helgeson  et  a 

.  1984 

0 

16 

84 

25 

7 

-8.65 

Helgeson  et  a 

.  1984 

0 

16 

84 

150 

8.94 

-7.79 

Gautier  et  al.  1994 

0 

16 

84 

150 

8.94 

-7.6 

Gautier  et  al.  1994 

~ 

~ 

97 

300 

9 

-6.3 

Alekseyev  et  al.  1997 

— 

— 

97 

300 

9 

-6.3 

Alekseyev  et  al.  1997 

25 


Table  14.  K-feldspar  Dissolution  rate  Data. 


an 

Composition 

ab 

or 

T,  °C 

bPH 

clog  k 

Source 

0 

7 

93 

70 

9 

-12.03 

Bevan  and  Savage  1989 

0 

7 

93 

95 

9 

-11.63 

Bevan  and  Savage  1989 

0 

7 

93 

95 

9 

-11.72 

Bevan  and  Savage  1989 

0 

6 

95 

70 

9.51 

-10.7 

Knauss  and  Copenhaver  1995 

a.  -  =  Not  given. 

b.  Acid  mechanism,  pH  4-5.91;  Neutral  mechanism,  pH  6-7;  base  mechanism,  pH  8.94-13. 

c.  Log  mole  m'2  s'1. 


Table  15.  K-Feldspar  Dissolution  Rate  Parameters. 


Acid  Mechanism 

Neutral  Mechanism 

Base  Mechanism 

alog  k  bE  cn 

alog  k 

bE 

alog  k 

bE  cn 

-10.06  51.7  0.500 

-12.41 

38.0 

-21.20 

94.1  -0.823 

a.  Rate  constant  k  computed  from  A  and  E,  25°C,  pH 

=  0,  mole  m'2  s'1. 

b.  Arrhenius  activation  energy  E,  kJ  mole 

c.  Reaction  order  n  with  respect  to  H+. 


3.1.3  Feldspathoids:  Nepheline  and  Leucite 

Dissolution  rate  data  for  nepheline  were  evaluated  from  Brady  and  Walther  (1989), 
Hamilton  et  al.  (2001),  Murphy  and  Helgeson  (1989),  Sverdrup  (1990),  and  Tole  et  al.  (1986), 
and  for  leucite,  from  Sverdrup  (1990).  Data  for  nepheline  were  selected  from  Tole  et  al.  (1993) 
for  the  acidic  and  neutral  mechanisms  because  they  were  obtained  from  experiments  at  a  wider 
range  of  temperatures  than  the  other  references.  The  data  from  Tole  et  al.  (1993)  are 
summarized  in  Table  16,  and  the  rate  parameters  determined  by  linear  regression  herein  using 
eqns.  (11)  and  (12)  are  given  in  Table  18.  Rate  parameters  for  the  basic  mechanism  for 
nepheline  and  acidic,  neutral,  and  basic  mechanisms  for  leucite  were  compiled  from  Sverdrup 
(1990)  and  are  summarized  in  Table  17,  with  activation  energies  calculated  herein  from  reported 
rate  constants  at  8  and  25  °C;  the  results  are  also  tabulated  in  Table  18. 
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Table  16.  Nepheline  Dissolution  Rate  Data. 


T,  °C 

apH 

blog  k 

Source 

25 

3 

-6.30 

Tole  et  al.  1986 

25 

3 

-6.30 

Tole  et  al.  1986 

60 

3 

-4.89 

Tole  et  al.  1986 

60 

3 

-4.89 

Tole  et  al.  1986 

80 

3 

-4.33 

Tole  et  al.  1986 

80 

3 

-4.33 

Tole  et  al.  1986 

25 

5 

-8.24 

Tole  et  al.  1986 

25 

5 

-8.24 

Tole  et  al.  1986 

60 

5 

-7.30 

Tole  et  al.  1986 

60 

5 

-7.30 

Tole  et  al.  1986 

80 

5 

-6.78 

Tole  et  al.  1986 

80 

5 

-6.78 

Tole  et  al.  1986 

25 

7 

-8.62 

Tole  et  al.  1986 

25 

7 

-8.62 

Tole  et  al.  1986 

60 

7 

-7.41 

Tole  et  al.  1986 

60 

7 

-7.41 

Tole  et  al.  1986 

80 

7 

-6.89 

Tole  et  al.  1986 

25 

11 

-8.56 

Tole  et  al.  1986 

60 

11 

-7.14 

Tole  et  al.  1986 

a.  Acid  mechanism,  pH  3-5;  Neutral  mechanism,  pH  7-11. 

b.  Log  mole  rrf2  s'1. 


Table  17.  Leucite  Dissolution  Rate  Constants  and  Reaction  Orders3. 


Acid  Mechanism 

Neutral  Mechanism 

Base  Mechanism 

log  k8»c  log  k25°c  bn 

log  k8°c  log  k25»c 

log  k8°c 

log  k25°c  bn 

Leucite  -10.40  -9.0  0.700 

-13.00  -12.2 

-11.40 

-10.8  -0.200 

a.  Data  from  Sverdrup,  1990;  log  k  units,  log  mole  m'2  s'1. 

b.  Reaction  order  n  with  respect  to  H+. 

Table  18.  Feldspathoid  Dissolution  Rate  Parameters. 

Acid  Mechanism 

Neutral  Mechanism 

Base  Mechanism 

alog  k  bE  cn 

alog  k  bE 

alog  k 

bE  cn 

Nepheline  -2.73  62.9  1.130  -8.56  65.4  -10.76  37.8  -0.200 

Leucite  -6.00  132.2  0.700  -9.20  75.5  -10.66  56.6  -0.200 


a.  Rate  constant  k  computed  from  A  and  E,  25°C,  pH  =  0,  mole  m'2  s'1. 

b.  Arrhenius  activation  energy  E,  kJ  mole'1. 

c.  Reaction  order  n  with  respect  to  H+. 
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3.2  Orthosilicates 


3.2.1  Olivine  Group 

Dissolution  rate  data  for  olivine  were  evaluated  from  Awad  et  al.  (2000),  Bailey  (1974), 
Blum  and  Lasaga  (1986),  Chen  and  Brantley  (2000),  Grandstaff  (1986),  Jonckbloedt  (1998), 

Luce  (1972),  Murphy  and  Helgeson  (1987;  1989),  Oelkers  (2001a),  Pokrovsky  and  Schott 
(2000b),  Pokrovsky  and  Schott  (2000a),  Rosso  and  Rimstidt  (2000),  Sverdrup  (1990),  Van  Herk 
et  al.  (1989),  Westrich  et  al.  (1993),  Wogelius  and  Walther,  (1991),  Wogelius  and  Walther 
(1992)  and  Xiao  et  al.  (1999). 

Dissolution  rate  data  were  selected  for  forsterite  from  Chen  and  Brantley  (2000),  Oelkers 
(2001a),  Pokrovsky  and  Schott  (2000b),  Pokrovsky  and  Schott  (2000a),  and  Wogelius  and 
Walther  (1992);  these  data  are  summarized  in  Table  19.  The  acidic  mechanism  for  forsterite 
apparently  operates  at  pH  up  to  ~9,  and  there  apparently  is  no  basic  mechanism.  Rate  parameters 
for  the  acidic  and  neutral  mechanisms  were  regressed  using  eqn.  (11),  and  the  results  are 
summarized  in  Table  20.  Rate  parameters  for  the  acidic  mechanism  for  fayalite  is  compiled  from 
Sverdrup  (1990;  Table  20),  with  activation  energies  calculated  herein  from  reported  rate 
constants  at  8  and  25  °C;  for  the  neutral  mechanism,  a  “cut-off’  rate  constant  at  25  °C  was 
calculated  so  that  rates  do  not  decrease  with  pH  increasing  above  9,  and  the  activation  energy  is 
set  equal  to  the  acidic  mechanism.  The  rate  parameters  for  fayalite  are  also  summarized  in  Table 
20. 


Table  19.  Forsterite  Dissolution  Rate  Data. 


Composition 

T,  °C 

apH 

blog  k 

Source 

fo91 

25 

1.03 

-7.23 

Pokrovsky  and  Schott,  2000 

fo91 

25 

1.11 

-7.33 

Pokrovsky  and  Schott,  2000 

fo91 

25 

1.12 

-7.28 

Pokrovsky  and  Schott,  2000 

fo91 

65 

1.8 

-6.29 

Wogelius  and  Walther,  1992 

fo91 

25 

2 

-7.80 

Oelkers,  2001 

fo91 

25 

2 

-7.75 

Oelkers,  2001 

fo91 

25 

2 

-7.92 

Oelkers,  2001 

fo91 

25 

2 

-7.91 

Oelkers,  2001 

fo91 

25 

2 

-7.85 

Oelkers,  2001 

fo91 

25 

2 

-7.81 

Oelkers,  2001 

fo91 

25 

2 

-7.98 

Oelkers,  2001 

fo91 

25 

2 

-7.88 

Oelkers,  2001 

fo91 

25 

2 

-7.88 

Oelkers,  2001 

fo91 

25 

2 

-7.90 

Oelkers,  2001 

fo91 

25 

2 

-7.89 

Oelkers,  2001 

fo91 

25 

2 

-7.92 

Oelkers,  2001 

fo91 

25 

2 

-7.90 

Oelkers,  2001 

fo91 

25 

2 

-7.91 

Oelkers,  2001 
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Table  19.  Forsterite  Dissolution  Rate  Data  -  continued. 


Composition 

T,  °C 

apH 

blog  k 

Source 

fo91 

25 

2 

-7.87 

Oelkers,  2001 

fo91 

25 

2 

-7.93 

Oelkers,  2001 

fo91 

25 

2 

-7.91 

Oelkers,  2001 

fo91 

25 

2 

-7.93 

Oelkers,  2001 

fo91 

25 

2 

-7.89 

Oelkers,  2001 

fo91 

25 

2 

-7.99 

Oelkers,  2001 

fo91 

25 

2 

-7.89 

Oelkers,  2001 

fo91 

25 

2 

-7.92 

Oelkers,  2001 

fo91 

35 

2 

-7.58 

Oelkers,  2001 

fo91 

45 

2 

-7.03 

Oelkers,  2001 

fo91 

55 

2 

-6.97 

Oelkers,  2001 

fo91 

65 

2 

-6.09 

Chen  and  Brantley,  2000 

fo91 

65 

2 

-6.57 

Oelkers,  2001 

fo91 

25 

2.13 

-7.70 

Pokrovsky  and  Schott,  2000 

fo91 

25 

2.18 

-7.71 

Pokrovsky  and  Schott,  2000 

fo91 

25 

2.21 

-7.72 

Pokrovsky  and  Schott,  2000 

fo91 

25 

2.7 

-7.73 

Pokrovsky  and  Schott,  2000 

fo91 

25 

2.74 

-7.83 

Pokrovsky  and  Schott,  2000 

fo91 

25 

2.81 

-7.82 

Pokrovsky  and  Schott,  2000 

fo91 

25 

2.82 

-7.89 

Pokrovsky  and  Schott,  2000 

fo91 

25 

2.85 

-8.30 

Pokrovsky  and  Schott,  2000 

f091 

65 

2.95 

-6.76 

Chen  and  Brantley,  2000 

f091 

65 

2.95 

-6.83 

Chen  and  Brantley,  2000 

f091 

65 

2.95 

-6.85 

Chen  and  Brantley,  2000 

f091 

65 

2.95 

-6.82 

Chen  and  Brantley,  2000 

fo91 

65 

2.95 

-6.68 

Chen  and  Brantley,  2000 

fo91 

25 

2.99 

-7.93 

Pokrovsky  and  Schott,  2000 

fo91 

25 

3.03 

-7.97 

Pokrovsky  and  Schott,  2000 

fo91 

25 

3.05 

-8.18 

Pokrovsky  and  Schott,  2000 

fo91 

25 

3.05 

-8.28 

Pokrovsky  and  Schott,  2000 

fo91 

25 

3.08 

-8.14 

Pokrovsky  and  Schott,  2000 

fo91 

25 

3.16 

-8.25 

Pokrovsky  and  Schott,  2000 

fo91 

25 

3.26 

-8.34 

Pokrovsky  and  Schott,  2000 

fo91 

25 

3.34 

-8.55 

Pokrovsky  and  Schott,  2000 

fo91 

25 

3.6 

-8.80 

Pokrovsky  and  Schott,  2000 

f091 

65 

4 

-7.60 

Chen  and  Brantley,  2000 

fo91 

25 

4.15 

-9.03 

Pokrovsky  and  Schott,  2000 

fo91 

25 

4.2 

-8.69 

Pokrovsky  and  Schott,  2000 

fo91 

25 

4.22 

-8.78 

Pokrovsky  and  Schott,  2000 

fo91 

25 

4.41 

-8.80 

Pokrovsky  and  Schott,  2000 

fo91 

25 

4.49 

-8.98 

Pokrovsky  and  Schott,  2000 
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Table  19.  Forsterite  Dissolution  Rate  Data  -  continued. 


Composition 

T,  °C 

apH 

blog  k 

Source 

fo91 

25 

4.55 

-9.03 

Pokrovsky  and  Schott,  2000 

fo91 

25 

4.77 

-8.98 

Pokrovsky  and  Schott,  2000 

fo91 

25 

4.85 

-9.02 

Pokrovsky  and  Schott,  2000 

fo91 

25 

4.95 

-9.15 

Pokrovsky  and  Schott,  2000 

f091 

65 

5 

-7.92 

Chen  and  Brantley,  2000 

fo91 

25 

5.4 

-9.31 

Pokrovsky  and  Schott,  2000 

fo91 

25 

5.44 

-9.47 

Pokrovsky  and  Schott,  2000 

fo91 

25 

5.44 

-9.47 

Pokrovsky  and  Schott,  2000 

fo91 

25 

5.47 

-9.57 

Pokrovsky  and  Schott,  2000 

fo91 

25 

5.61 

-9.79 

Pokrovsky  and  Schott,  2000 

fo91 

25 

5.7 

-9.49 

Pokrovsky  and  Schott,  2000 

fo91 

65 

6 

-8.13 

Wogelius/Walther,  1992 

fo91 

25 

6.13 

-9.92 

Pokrovsky  and  Schott,  2000 

fo91 

25 

6.18 

-9.60 

Pokrovsky  and  Schott,  2000 

fo91 

25 

6.28 

-9.72 

Pokrovsky  and  Schott,  2000 

fo91 

25 

6.39 

-10.04 

Pokrovsky  and  Schott,  2000 

fo91 

25 

6.6 

-9.97 

Pokrovsky  and  Schott,  2000 

fo91 

25 

7.19 

-10.11 

Pokrovsky  and  Schott,  2000 

fo91 

25 

7.21 

-10.04 

Pokrovsky  and  Schott,  2000 

fo91 

25 

7.25 

-9.88 

Pokrovsky  and  Schott,  2000 

fo91 

25 

7.33 

-9.93 

Pokrovsky  and  Schott,  2000 

fo91 

25 

7.6 

-10.47 

Pokrovsky  and  Schott,  2000 

fo91 

25 

7.9 

-10.57 

Pokrovsky  and  Schott,  2000 

fo91 

25 

7.98 

-10.39 

Pokrovsky  and  Schott,  2000 

fo91 

25 

8.02 

-10.13 

Pokrovsky  and  Schott,  2000 

fo91 

25 

8.13 

-10.33 

Pokrovsky  and  Schott,  2000 

fo91 

25 

8.38 

-10.18 

Pokrovsky  and  Schott,  2000 

fo91 

25 

8.5 

-10.33 

Pokrovsky  and  Schott,  2000 

fo91 

25 

8.53 

-10.17 

Pokrovsky  and  Schott,  2000 

fo91 

25 

8.55 

-10.42 

Pokrovsky  and  Schott,  2000 

fo91 

25 

8.55 

-10.40 

Pokrovsky  and  Schott,  2000 

fo91 

25 

8.63 

-10.44 

Pokrovsky  and  Schott,  2000 

fo91 

25 

8.68 

-10.46 

Pokrovsky  and  Schott,  2000 

fo91 

25 

8.71 

-10.60 

Pokrovsky  and  Schott,  2000 

fo91 

25 

9.3 

-10.49 

Pokrovsky  and  Schott,  2000 

fo91 

25 

9.58 

-10.54 

Pokrovsky  and  Schott,  2000 

fo91 

25 

10.08 

-10.79 

Pokrovsky  and  Schott,  2000 

fo91 

25 

10.24 

-10.72 

Pokrovsky  and  Schott,  2000 

fo91 

25 

10.76 

-10.38 

Pokrovsky  and  Schott,  2000 

fo91 

25 

10.78 

-10.95 

Pokrovsky  and  Schott,  2000 

fo91 

25 

10.85 

-10.31 

Pokrovsky  and  Schott,  2000 
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Table  19.  Forsterite  Dissolution  Rate  Data  -  continued. 


Composition 

T,  °C 

apH 

blog  k 

Source 

fo91 

25 

10.9 

-10.44 

Pokrovsky  and  Schott,  2000 

fo91 

25 

10.9 

-10.69 

Pokrovsky  and  Schott,  2000 

fo91 

25 

10.95 

-10.60 

Pokrovsky  and  Schott,  2000 

fo91 

25 

11.01 

-10.57 

Pokrovsky  and  Schott,  2000 

fo91 

25 

11.1 

-10.65 

Pokrovsky  and  Schott,  2000 

fo91 

25 

11.11 

-10.56 

Pokrovsky  and  Schott,  2000 

fo91 

25 

11.12 

-10.52 

Pokrovsky  and  Schott,  2000 

fo91 

25 

11.13 

-10.71 

Pokrovsky  and  Schott,  2000 

fo91 

25 

11.15 

-10.57 

Pokrovsky  and  Schott,  2000 

fo91 

25 

12.06 

-10.85 

Pokrovsky  and  Schott,  2000 

a.  Acid  mechanism,  pH  1.03-8.71;  Neutral  mechanism,  pH  9.3-12.06. 

b.  Log  mole  m"2  s'1. 


Table  20.  Orthosilicate  Dissolution  Rate  Constants  and  Reaction  Orders3. 


Acid  Mechanism 

Neutral  Mechanism 

Base  Mechanism 

log  k8»c 

log  k25°c 

bn 

log  k8°c 

log  k25"C 

log  k8°c 

log  k25°c  bn 

Fayalite 

-5.8 

-4.8 

1.0 

- 

- 

- 

- 

Almandine 

-6.2 

-5.2 

1.0 

-11.8 

-10.7 

-9.1 

-8.7  -0.3 

Grossular 

-6 

-5.1 

1.0 

-11.8 

-10.7 

— 

.. 

Andradite 

-6.2 

-5.2 

1.0 

-11.8 

-10.7 

— 

.. 

Staurolite 

-7.1 

-6.9 

1.0 

-12.8 

-12.2 

-11.1 

-10.6  -0.3 

Zoisite 

-8.2 

-7.5 

0.5 

- 

- 

~ 

~ 

a.  Data  from  Sverdrup,  1990;  log  k  units,  log  mole  m'2  s'1. 

b.  Reaction  order  n  with  respect  to  H+. 


3.2.2  Garnet  Group 

Rate  parameters  for  the  acidic,  neutral,  and  basic  mechanisms  for  almandine,  and  acidic 
and  neutral  mechanisms  (data  are  absent  for  the  basic  mechanism)  for  grossular  and  andradite  are 
compiled  from  Sverdrup  (1990;  Table  20),  with  activation  energies  calculated  herein  from 
reported  rate  constants  at  8  and  25  °C.  The  rate  parameters  for  the  garnet  group  minerals  are 
summarized  in  Table  22. 
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3.2.3  AI2Si05  Group 


Dissolution  rate  data  for  kyanite  were  obtained  from  Oelkers  and  Schott  (1999;  Table  21) 
for  the  acidic  mechanism  only.  The  rate  parameters  are  regressed  herein  using  eqn.  (11),  and  are 
summarized  in  Table  23.  For  the  neutral  mechanism,  a  “cut-off’  rate  constant  at  25  °C  was 
calculated  so  that  rates  do  not  decrease  with  pH  increasing  above  7,  and  the  activation  energy  is 
set  equal  to  that  for  the  acidic  mechanism. 

Table  21 .  Kyanite  Dissolution  Rate  Data. 


T,  °C 

pH 

alog  k 

T,  °C 

pH 

alog  k 

150 

1.65 

-9.63 

150 

2.04 

-9.87 

150 

1.65 

-9.55 

150 

2.04 

-9.87 

150 

1.65 

-9.42 

150 

2.04 

-9.86 

150 

1.85 

-9.84 

150 

2.04 

-9.79 

150 

1.85 

-9.66 

150 

2.04 

-9.88 

150 

1.85 

-9.65 

150 

2.04 

-9.79 

108 

2.03 

-10.96 

150 

2.04 

-9.82 

108 

2.03 

-10.96 

131 

2.04 

-9.98 

111 

2.03 

-10.54 

174 

2.06 

-9.89 

113 

2.03 

-10.28 

175 

2.06 

-9.67 

121 

2.03 

-10.48 

175 

2.06 

-9.69 

121 

2.03 

-10.48 

175 

2.06 

-9.62 

121 

2.03 

-10.66 

175 

2.06 

-9.69 

122 

2.03 

-10.48 

175 

2.06 

-9.61 

125 

2.03 

-10.40 

175 

2.06 

-9.56 

125 

2.03 

-10.34 

175 

2.06 

-9.77 

125 

2.03 

-10.28 

194 

2.07 

-9.46 

126 

2.03 

-10.40 

194 

2.07 

-9.27 

126 

2.03 

-10.29 

150 

2.24 

-10.40 

150 

2.04 

-9.86 

150 

2.24 

-10.24 

150 

2.04 

-9.92 

150 

2.24 

-10.21 

150 

2.04 

-9.84 

150 

2.24 

-10.37 

a.  Log  mole  rrf2  s'1. 
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3.2.4  Staurolite 


Rate  parameters  for  the  acidic, 
compiled  from  Sverdrup  (1990;  Table 
reported  rate  constants  at  8  and  25  °C. 
Table  23. 


neutral,  and  basic  mechanisms  for  staurolite  were 
20),  with  activation  energies  calculated  herein  from 
The  rate  parameters  for  the  staurolite  are  summarized  in 


3.2.5  Epidote  Group 

Dissolution  rate  data  for  epidote  and  zoisite  were  evaluated  from  Rose  (1991), 
Kalinowski  et  al.  (1998),  and  Sverdrup  (1990).  Data  for  the  acidic,  neutral,  and  basic 
mechanisms  for  epidote  from  Rose  (1991),  and  Kalinowski  et  al.  (1998)  were  selected,  and  are 
tabulated  in  Table  22.  The  rate  parameters  for  epidote  are  regressed  herein  using  eqn.  (11),  and 
are  summarized  in  Table  23.  Rate  parameters  for  the  acidic  mechanism  for  zoisite  was  compiled 
from  Sverdrup  (1990;  Table  20),  with  the  activation  energy  calculated  herein  from  reported  rate 
constants  at  8  and  25  °C;  for  the  neutral  mechanism,  a  “cut-off’  rate  constant  at  25  °C  was 
calculated  so  that  rates  do  not  decrease  with  pH  increasing  above  7,  and  the  activation  energy  is 
set  equal  to  the  acidic  mechanism.  The  rate  parameters  for  zoisite  are  summarized  in  Table  23. 
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Table  22.  Epidote  Dissolution  Rate  Data. 


T,  °C 

apH 

blog  k 

Source 

25 

0.2 

-10.23 

Kalinowski  et  al.  1998 

25 

1.08 

-10.49 

Kalinowski  et  al.  1998 

25 

1.29 

-11.31 

Kalinowski  et  al.  1998 

90 

1.29 

-8.715 

Rose  1991 

25 

1.4 

-11.53 

Kalinowski  et  al.  1998 

25 

1.4 

-11.528 

Rose  1991 

50 

1.4 

-10.322 

Rose  1991 

70 

1.4 

-9.268 

Rose  1991 

70 

1.4 

-9.63 

Rose  1991 

70 

1.4 

-9.677 

Rose  1991 

90 

1.4 

-8.9 

Rose  1991 

25 

1.99 

-11.19 

Kalinowski  et  al.  1998 

25 

2.01 

-10.73 

Kalinowski  et  al.  1998 

25 

2.5 

-11.8 

Kalinowski  et  al.  1998 

90 

2.5 

-9.2 

Rose  1991 

90 

2.5 

-9.269 

Rose  1991 

25 

2.98 

-11.36 

Kalinowski  et  al.  1998 

25 

2.99 

-11.63 

Kalinowski  et  al.  1998 

25 

3 

-11.64 

Kalinowski  et  al.  1998 

25 

3.03 

-11.42 

Kalinowski  et  al.  1998 

25 

3.47 

-12.02 

Kalinowski  et  al.  1998 

90 

3.47 

-9.424 

Rose  1991 

90 

3.47 

-9.424 

Rose  1991 

90 

3.47 

-9.528 

Rose  1991 

25 

4.02 

-12.19 

Kalinowski  et  al.  1998 

25 

4.05 

-11.72 

Kalinowski  et  al.  1998 

25 

4.05 

-12 

Kalinowski  et  al.  1998 

25 

4.12 

-11.65 

Kalinowski  et  al.  1998 

25 

4.38 

-12.17 

Kalinowski  et  al.  1998 

90 

4.38 

-9.578 

Rose  1991 

25 

5.5 

-12.15 

Kalinowski  et  al.  1998 

90 

5.5 

-9.55 

Rose  1991 

25 

5.6 

-11.8 

Kalinowski  et  al.  1998 

25 

6.48 

-12.29 

Kalinowski  et  al.  1998 

25 

6.5 

-11.93 

Kalinowski  et  al.  1998 

90 

6.62 

-9.795 

Rose  1991 

90 

7.75 

-9.892 

Rose  1991 

90 

8.64 

-10.048 

Rose  1991 

90 

10.83 

-8.831 

Rose  1991 

25 

10.6 

-11.44 

Kalinowski  et  al.  1998 

a.  Acid  mechanism,  pH  1.2  -4.38;  Neutral  mechanism,  pH  5.5-7.75;  basic 
mechanism,  8.64-10.8. 

b.  Log  mole  m'2  s'1. 
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Table  23.  Orthosilicate  Dissolution  Rate  Parameters. 


Acid  Mechanism 

Neutral  Mechanism 

Base  Mechanism 

alog  k 

bE 

cn 

alog  k 

bE 

alog  k 

bE 

cn 

Forsterite 

-6.85 

67.2 

0.470 

-10.64 

79.0 

— 

— 

— 

Fayalite 

-4.80 

94.4 

-12.80 

94.4 

- 

- 

- 

Almandine 

-5.20 

94.4 

1.000 

-10.70 

103.8 

-13.71 

37.8 

-0.350 

Grossular 

-5.10 

85.0 

1.000 

-10.70 

103.8 

— 

— 

— 

Andradite 

-5.20 

94.4 

1.000 

-10.70 

103.8 

— 

— 

— 

Kyanite 

-10.17 

-53.9 

1.268 

-17.44 

53.9 

~ 

~ 

~ 

Staurolite 

-6.90 

18.9 

1.000 

-12.20 

56.6 

-14.90 

47.2 

-0.300 

Epidote 

-10.60 

71.1 

0.338 

-11.99 

70.7 

-17.33 

79.1 

-0.556 

Zoisite 

-7.50 

66.1 

0.500 

-11.20 

66.1 

~ 

~ 

~ 

a.  Rate  constant  k  computed  from  A  and  E,  25°C,  pH  =  0,  mole  m'2  s'1. 

b.  Arrhenius  activation  energy  E,  kJ  mole'1. 

c.  Reaction  order  n  with  respect  to  H+. 


3.3  Cyclosilicates 

3.3.1  Cordierite  and  Tourmaline 

Rate  parameters  for  the  acidic  and  neutral  mechanisms  (data  are  absent  for  the  basic 
mechanism)  for  cordierite  and  tourmaline  were  compiled  from  Sverdrup  (1990;  Table  24),  with 
activation  energies  calculated  herein  from  reported  rate  constants  at  8  and  25  °C.  The  rate 
parameters  for  the  garnet  group  minerals  are  summarized  in  Table  22. 


Table  24.  Cyclosilicate  Dissolution  Rate  Constants  and  Reaction  Orders3. 


Acid  Mechanism 

Neutral  Mechanism 

log  k8°c  log  k25»c  bn 

log  k8»c 

log  k25°c 

Cordierite 

-5.0  -3.8  1.000 

-11.5 

-11.2 

Tourmaline 

-7.3  -6.5  1.000 

-12.1 

-11.2 

a.  Data  from  Sverdrup,  1990;  log  k  units,  log  mole  m'2  s'1. 

b.  Reaction  order  n  with  respect  to  H+. 
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Table  25.  Cyclosilicate  Dissolution  Rate  Parameters. 


Acid  Mechanism 

Neutral  Mechanism 

alog  k 

bE 

cn 

alog  k 

bE 

Cordierite 

-3.80 

113.3 

1.000 

-11.20 

28.3 

Tourmaline 

-6.50 

75.5 

1.000 

-11.20 

85.0 

a.  Rate  constant  k  computed  from  A  and  E,  25°C,  pH  =  0,  mole  m'2  s'1. 

b.  Arrhenius  activation  energy  E,  kJ  mole'1. 

c.  Reaction  order  n  with  respect  to  H+. 


3.4  Inosilicates 

3.4.1  Pyroxene  Group  /  Pyroxenoid  Group 

Dissolution  rate  data  for  pyroxenes  and  pyroxenoids  were  evaluated  from  Bailey  (1974), 
Brantley  and  Chen  (1995),  Chen  and  Brantley  (1998),  Knauss  et  al.  (1993),  Murphy  and 
Helgeson  (1987;  1989),  Rimstidt  and  Dove  (1986),  Schott  and  Berner  (1985),  Schott  et  al. 

(1981),  Siegel  and  Pfannkuch  (1984),  Sverdrup  (1990),  and  Weissbart  (2000). 

Dissolution  rate  data  for  the  acidic  and  neutral  mechanisms  (data  are  absent  for  the  basic 
mechanism)  were  selected  as  follows:  augite,  from  Schott  et  al.  (1981)  and  Sverdrup  (1990); 
bronzite,  from  Sverdrup  (1990);  diopside,  from  Brantley  and  Chen  (1995),  Chen  and  Brantley 
(1998),  Knauss  et  al.  (1993),  and  Schott  and  Berner  (1985);  enstatite,  from  Schott  and  Bemer 
(1985);  jadeite,  from  Sverdrup  (1990);  spodumene,  from  Sverdrup  (1990);  wollastonite,  from 
Murphy  and  Helgeson  (1987;  1989)  and  Rimstidt  and  Dove  (1986).  The  rate  parameters  were 
regressed  using  eqn.  (11),  except  for  the  bronzite  data  from  Sverdrup  (1990),  for  which  rate 
constants  at  8  and  25  °C  for  single  mechanisms  are  used  herein  to  calculate  the  activation  energy. 
The  resulting  rate  parameters  are  summarized  in  Table  26. 

For  enstatite  and  augite,  the  available  data  are  limited  to  pH  less  than  ~7.  Although  rates 
may  decrease  further  with  further  pH  increase  above  7,  a  “cut-off’  rate  constant  at  25  °C  was 
calculated  for  the  neutral  mechanism  so  that  rates  do  not  decrease  with  pH  increasing  above  7, 
and  the  activation  energy  is  set  equal  to  the  acidic  mechanism. 
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Table  26.  Pyroxene  and  Pyroxenoid  Dissolution  Rate  Parameters. 


Acid  Mechanism 

Neutral  Mechanism 

alog  k 

bE 

cn 

alog  k 

bE 

augite 

-6.82 

78.0 

0.700 

-11.97 

78.0 

bronzite 

-8.30 

47.2 

0.650 

-11.70 

66.1 

diopside 

-6.36 

96.1 

0.710 

-11.11 

40.6 

enstatite 

-9.02 

80.0 

0.600 

-12.72 

80.0 

jadeite 

-6.00 

132.2 

0.700 

-9.50 

94.4 

spodumene 

-4.60 

94.4 

0.700 

-9.30 

66.1 

wollastonite 

-5.37 

54.7 

0.400 

-8.88 

54.7 

a.  Rate  constant  k  computed  from  A  and  E,  25°C,  pH  =  0,  mole  m'2  s'1. 

b.  Arrhenius  activation  energy  E,  kJ  mole'1. 

c.  Reaction  order  n  with  respect  to  H+. 


3.4.2  Amphibole  Group 


Dissolution  rate  data  for  amphiboles  were  evaluated  from  Brantley  and  Chen  (1995), 
Chen  and  Brantley  (1998),  Sverdrup  (1990),  and  Zhang  and  Bloom  (1999).  Data  for  the  acidic 
and  neutral  mechanisms  (data  are  absent  for  the  basic  mechanism)  were  selected  as  follows: 
anthophyllite,  from  Chen  and  Brantley  (1998);  glaucophane,  hornblende,  riebeckite  and 
tremolite,  from  Sverdrup  (1990).  The  rate  parameters  were  regressed  using  eqn.  (11),  except  for 
the  data  from  Sverdrup  (1990),  for  which  rate  constants  at  8  and  25  °C  for  single  mechanisms  are 
used  herein  to  calculate  the  activation  energy.  The  resulting  rate  parameters  are  summarized  in 
Table  27. 

Table  27.  Amphibole  Dissolution  Rate  Parameters. 


Acid  Mechanism 

Neutral  Mechanism 

alog  k 

bE 

cn 

alog  k 

bE 

anthophyllite 

-11.94 

51.0 

0.440 

-14.24 

51.0 

glaucophane 

-5.60 

85.0 

0.700 

-10.10 

94.4 

hornblende 

-7.00 

75.5 

0.600 

-10.30 

94.4 

riebeckite 

-7.70 

56.6 

0.700 

-12.20 

47.2 

tremolite 

-8.40 

18.9 

0.700 

-10.60 

94.4 

a.  Rate  constant  k  computed  from  A  and  E,  25°C,  pH  =  0,  mole  m'2  s'1. 

b.  Arrhenius  activation  energy  E,  kJ  mole'1. 

c.  Reaction  order  n  with  respect  to  H+. 
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3.5  Phyllosilicates 


3.5.1  Mica  Group 

Dissolution  rate  data  for  the  micas  were  evaluated  from  Acker  and  Bricker  (1992), 
Clemency  and  Lin  (1981),  Kalinowski  and  Schweda  (1996),  Knauss  and  Wolery  (1989), 
Kuwahara  and  Aoke  (1999),  Lin  and  Clemency  (1981a;  1981b),  Malmstrom  and  Banwart 
(1997),  Nagy  (1995),  Robie  (1976),  Sverdrup  (1990),  Swoboda-Colberg  and  Drever  (1993), 
Taylor  et  al.  (2000b),  and  Trotignon  and  Turpault  (1992). 

Data  were  selected  as  follows:  biotite,  from  Acker  and  Bricker  (1992)  and  Nagy  (1995); 
glauconite,  from  Sverdrup  (1990);  muscovite,  paragonite,  phlogopite  and  pyrophyllite,  from 
Nagy  (1995).  The  rate  parameters  were  regressed  using  eqn.  (11),  except  for  the  data  from 
Sverdrup  (1990),  for  which  rate  constants  at  8  and  25  °C  for  single  mechanisms  are  used  herein 
to  calculate  the  activation  energy.  The  compiled  rate  parameters  are  summarized  in  Table  28. 


Table  28.  Orthosilicate  Dissolution  Rate  Parameters. 


Acid  Mechanism 

Neutral  Mechanism 

Base  Mechanism 

alog  k 

bE 

cn 

alog  k 

bE 

alog  k 

bE  cn 

biotite 

-9.84 

22.0 

0.525 

-12.55 

22.0 

— 

.. 

glauconite 

-4.80 

85.0 

0.700 

-9.10 

85.0 

- 

- 

muscovite 

-11.85 

22.0 

0.370 

-13.55 

22.0 

-14.55 

22.0  -0.220 

muscovite 

— 

- 

— 

-13.00 

22.0 

— 

— 

paragonite 

- 

- 

- 

-13.00 

22.0 

- 

- 

phlogopite 

- 

- 

- 

-12.40 

29.0 

- 

-- 

pyrophyllite 

- 

-- 

- 

-12.40 

29.0 

- 

-- 

a.  Rate  constant  k  computed  from  A  and  E,  25°C,  pH  =  0,  mole  m'2  s'1. 

b.  Arrhenius  activation  energy  E,  kJ  mole'1. 

c.  Reaction  order  n  with  respect  to  H+. 


3.5.2  Clay  Group 


Dissolution  rate  data  for  three  additional  clay  group  minerals  were  evaluated:  kaolinite 
(Bauer  and  Berger,  1998;  Cama  et  al.,  2002;  Carroll  and  Walther,  1990;  Carroll-Webb  and 
Walther,  1988;  Devidal  et  al.,  1997;  Ganor  et  al.,  1995;  Hayashi  and  Yamada,  1990;  Huertas  et 
al.,  1999a;  Huertas  et  al.,  2001b;  Huertas  et  al.,  1999b;  Metz  and  Ganor,  2001;  Nagy,  1995;  Nagy 
et  al.,  1991;  Nagy  and  Lasaga,  1990;  Nagy  and  Lasaga,  1993;  Oelkers  et  al.,  2001;  Polzer  and 
Hem,  1965;  Soong,  1993;  Soong  and  Bames,  1991;  Sutheimer  et  al.,  1999;  Wieland  and  Stumm, 
1992);  montmorillonite  (Furrer  et  al.,  1993;  Hayashi  and  Yamada,  1990;  Huertas  et  al.,  2001a; 
Nagy,  1995;  Zysset  and  Schindler,  1996);  and  smectite  (Altaner,  1986;  Bauer  and  Berger,  1998; 
Cama  et  al.,  2000;  Cama  et  al.,  1994;  Novak  and  Cicel,  1978). 

Data  for  the  clay  group  minerals  were  selected  as  follows:  kaolinite,  from  Carroll  and 
Walther  (1990),  Ganor  et  al.  (1995),  Huertas  et  al.  (1999a),  Huertas  et  al.  (1999b),  Nagy  et  al. 
(1991),  Soong  (1993);  montmorillonite,  from  Nagy  (1995);  smectite,  from  Bauer  and  Berger 
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(1998),  Huertas  et  al.  (2001a),  Sverdrup  (1990),  and  Zysset  and  Schindler  (1996).  The  rate 
parameters  were  regressed  using  eqn.  (11),  except  for  the  data  from  Sverdrup  (1990),  for  which 
rate  constants  at  8  and  25  °C  for  single  mechanisms  are  used  herein  to  calculate  the  activation 
energy.  The  compiled  rate  parameters  are  summarized  in  Table  29.  Nagy  et  al.  (1991)  have 
quantified  the  chemical  affinity  parameters  for  kaolinite  at  80  °C  and  pH  =  3;  with  the 
experimental  uncertainty,  the  experimental  data  can  be  fit  to  eqn.  (4)  with  either  p  =  1 .000  and  q 
=  0.897,  or p  =  0.850  and  q  =  1.000. 


Table  29.  Clay  Group  Mineral  Dissolution  Rate  Parameters. 


Acid  Mechanism 

Neutral  Mechanism 

Base  Mechanism 

alogk 

bE 

cn 

alog  k 

bE 

alog  k 

bE 

cn 

kaolinite  -11.31 

65.9 

0.777 

-13.18 

22.2 

-17.05 

17.9 

-0.472 

dmontmorillonite  -12.71 

48.0 

0.220 

-14.41 

48.0 

-14.41 

48. 

-0.130 

esmecite  -10.98 

23.6 

0.340 

-12.78 

35.0 

-16.52 

58.9 

-0.400 

a.  Rate  constant  k  computed  from  A  and  E,  25°C,  pH  =  0,  mole  m'2  s'1. 

b.  Arrhenius  activation  energy  E,  kJ  mole'1. 

c.  Reaction  order  n  with  respect  to  H+. 

d.  Montmorillonite,  Ko.3i8(Si3.975Alo.o25)(Ali.509Feo.205Mgo.283)(OH)2. 

e.  Smectite,  Ko.o4Cao.5(Al2.8Feo.53Mgo.7)(Si7.65Alo.35)02o(OH)4. 


3.5.3  Miscellaneous  Phyllosilicates 

Dissolution  rate  data  for  several  miscellaneous  phyllosilicate  minerals  were  evaluated: 
serpentine  (Bales  and  Morgan,  1985;  Cetisli  and  Gedikbey,  1990;  Gronow,  1987;  Hume  and 
Rimstidt,  1992;  Lin  and  Clemency,  1981b;  Nagy,  1995;  Sverdrup,  1990;  Thomassin  et  al.,  1977), 
chlorite  (Hayashi  and  Yamada,  1990;  Johnson  et  al.,  2001;  Nagy,  1995;  Sverdrup,  1990),  talc 
(Lin  and  Clemency,  1981b;  Nagy,  1995),  and  prehnite  (Rose,  1991). 

Data  for  these  phyllosilicate  minerals  were  selected  as  follows:  lizardite  serpentine,  from 
Sverdrup  (1990);  chrysotile  serpentine,  from  Nagy  (1995);  chlorite,  from  Nagy  (1995);  talc,  from 
Nagy  (1995);  and  prehnite,  from  Rose  (1991).  The  compiled  parameters  are  summarized  in 
Table  30. 
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Table  30.  Miscellaneous  Phyllosilicate  Mineral  Dissolution  Rate  Parameters. 


Acid  Mechanism 

Neutral  Mechanism 

Base  Mechanism 

alog  k 

bE 

cn 

alog  k 

bE 

alog  k 

bE 

cn 

lizardite 

-5.70 

75.5 

0.800 

-12.40 

56.6 

— 

— 

— 

chrysotile 

- 

- 

~ 

-12.00 

73.5 

-13.58 

73.5 

-0.230 

dchlorite 

-11.11 

88.0 

0.500 

-12.52 

88.0 

— 

— 

— 

talc 

- 

— 

— 

-12.00 

42.0 

— 

— 

— 

prehnite 

-10.66 

80.5 

0.256 

-13.16 

93.4 

-14.86 

93.4 

-0.200 

a.  Rate  constant  k  computed  from  A  and  E,  25°C,  pH  =  0,  mole  m'2  s'1. 

b.  Arrhenius  activation  energy  E,  kJ  mole'1. 

c.  Reaction  order  n  with  respect  to  H+. 

d.  Clinochlore  14  A. 


3.6  Oxides 


Dissolution  rate  data  for  several  oxide  minerals  were  evaluated:  goethite  (Ruan  and 
Gilkes,  1995),  hematite  (Bruno  and  Duro,  2000;  Bruno  et  al.,  1992;  Hummel,  2000;  Petrie,  1995; 
Ruan  and  Gilkes,  1995),  ilmenite  (White  et  al.,  1994),  magnetite  (White  et  al.,  1994),  manganite 
(Petrie,  1995),  pyrolusite  (Petrie,  1995),  and  uraninite  (Grandstaff,  1976).  Data  for  the  oxide 
minerals  were  selected  as  follows:  goethite,  from  Ruan  and  Gilkes  (1995);  hematite,  Ruan  and 
Gilkes  (1995);  ilmenite,  White  et  al.  (1994);  magnetite,  White  et  al.  (1994);  pyrolusite,  Petrie 
(1995);  and  uraninite  Grandstaff  (1976).  The  compiled  parameters  are  summarized  in  Table  31. 


Table  31.  Oxide  Mineral  Dissolution  Rate  Parameters. 


Acid  Mechanism 

Neutral  Mechanism 

alog  k 

bE 

cn 

alog  k 

bE 

goethite 

~ 

~ 

~ 

-7.94 

86.5 

hematite 

-9.39 

66.2 

1.000 

-14.60 

66.2 

magnetite 

-8.59 

18.6 

0.279 

-10.78 

18.6 

ilmenite 

-8.35 

37.9 

0.421 

-11.16 

37.9 

uraninite 

~ 

~ 

~ 

-7.98 

32.0 

a.  Rate  constant  k  computed  from  A  and  E,  25°C,  pH  =  0,  mole  m'2  s'1. 

b.  Arrhenius  activation  energy  E,  kJ  mole'1. 

c.  Reaction  order  n  with  respect  to  H+. 
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3.7  Hydroxides 


Dissolution  rate  data  for  brucite  (Jordan  and  Rammensee,  1996;  Lin  and  Clemency, 
1981b;  Nagy,  1995;  Sverdrup,  1990),  gibbsite  (Bloom,  1983;  Bloom  and  Erich,  1987;  Chang  et 
ah,  1979;  Ganor  et  ah,  1999;  Mogollon  et  ah,  1996;  Mogollon  et  ah,  2000;  Nagy,  1995;  Nagy 
and  Lasaga,  1990;  Nagy  and  Lasaga,  1992;  Nagy  and  Lasaga,  1993;  Packter  and  Dhillon,  1974; 
Scotford  and  Glastonbury,  1972a;  Scotford  and  Glastonbury,  1972b),  and  diaspore  (Chang  et  ah, 
1979)  were  evaluated.  Data  for  brucite  were  selected  from  Nagy  (1995)  and  Jordan  and 
Rammensee  (1996).  Data  for  gibbsite  were  selected  from  Nagy  and  Lasaga  (1992)  and  Packter 
and  Dhillon  (1974).  The  compiled  parameters  are  summarized  in  Table  32. 


Table  32.  Hydroxide  Mineral  Dissolution  Rate  Parameters. 


Acid  Mechanism 

Neutral  Mechanism 

Base  Mechanism 

alog  k 

bE 

cn 

alog  k 

bE 

alog  k 

bE 

cn 

brucite 

-4.73 

59.0 

0.500 

-8.24 

42.0 

— 

— 

— 

gibbsite 

-7.65 

47.5 

0.992 

-11.50 

61.2 

-16.65 

80.1 

-0.784 

diaspore 

~ 

~ 

~ 

-13.33 

47.5 

-23.60 

47.5 

-1.503 

a.  Rate  constant  k  computed  from  A  and  E,  25°C,  pH  =  0,  mole  m'2  s'1. 

b.  Arrhenius  activation  energy  E,  kJ  mole'1. 

c.  Reaction  order  n  with  respect  to  H+. 


3.8  Carbonates 


Dissolution  rate  data  for  four  carbonate  minerals  were  evaluated:  calcite  (Alkattan  et  ah, 
1998;  Alkattan  et  ah,  2002;  Bemer  and  Morse,  1974;  Brown  et  ah,  1993;  Busenberg  and 
Plummer,  1986;  Chou  et  al.,  1989;  Deleuze  and  Brantley,  1997;  House,  1981;  Inskeep  and 
Bloom,  1985;  Jeschke  and  Dreybrodt,  2002b;  Jimenez-Lopez  et  al.,  2001;  Jordan  and 
Rammensee,  1998;  Lebron  and  Suarez,  1996;  Lebron  and  Suarez,  1998;  Maclnnis  and  Brantley, 
1992;  Maldonado  et  al.,  1992;  Morse,  1978;  Mucci,  1986;  Mucci  and  Morse,  1983;  Plummer  and 
Wigley,  1976;  Plummer  et  al.,  1978;  Plummer  et  al.,  1979;  Reddy,  1986;  Reddy  and  Gaillard, 
1981;  Reddy  et  al.,  1981;  Schott  et  al.,  1989;  Shiraki  and  Brantley,  1995;  Sjoberg,  1976;  Sjoberg 
and  Rickard,  1983;  Sjoberg  and  Rickard,  1984a;  Sjoberg  and  Rickard,  1984b;  Sjoberg  and 
Rickard,  1985;  Talman  et  al.,  1990;  Teng  et  al.,  2000;  Zhang  and  Dawe,  1998;  Zhang  and  Dawe, 
2000;  Zhong  and  Mucci,  1989;  Zhong  and  Mucci,  1993;  Zuddas  and  Mucci,  1994;  Zuddas  and 
Mucci,  1998);  dawsonite  (Johnson  et  al.,  2001);  dolomite  (Arvidson  and  Mackenzie,  1999;  Baker 
and  Kastner,  1981a;  Baker  and  Kastner,  1981b;  Busenberg  and  Plummer,  1982;  Busenberg  and 
Plummer,  1986;  Chou  et  al.,  1989;  Gaines,  1980;  Gautelier  et  al.,  1999;  Katz  and  Matthews, 

1977;  Liittge  et  al.,  2003;  Nordeng  and  Sibley,  1994;  Sibley  et  al.,  1987);  magnesite  (Chou  et  al., 
1989;  Higgins  et  al.,  2002;  Johnson  et  al.,  2001;  Jordan  et  al.,  2001;  Pokrovsky  and  Schott, 

1999);  and  siderite  (Greenberg  and  Tomson,  1992;  Jensen  et  al.,  2002). 

Data  for  the  carbonate  minerals  were  selected  as  follows:  calcite,  from  Plummer  et  al. 
(1978)  and  Talman  et  al.  (1990);  dawsonite,  Johnson  et  al.  (2001);  dolomite,  Busenberg  and 
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Plummer  (1982);  and  magnesite,  Chou  et  al.  (1989)  and  Pokrovsky  and  Schott  (1999).  The  rate 
parameters  were  regressed  using  eqn.  (11),  and  are  summarized  in  Table  33.  In  the  absence  of 
data  for  the  activation  energy  for  magnesite  dissolution,  the  activation  energy  of  calcite  is 
assigned  to  magnesite.  Pokrovsky  and  Schott  (1999)  have  quantified  the  chemical  affinity 
parameters  for  magnesite  at  25  °C  with  pH  ranging  from  0.2  to  12,  with p  =  4.00  and  q  =  1.00. 


Table  33.  Carbonate  Mineral  Dissolution  Rate  Parameters. 


Acid  Mechanism 

Neutral  Mechanism 

Carbonate  Mechanism 

alog  k 

bE 

cn 

alog  k 

bE 

alog  k 

bE 

dn 

calcite 

-0.30 

14.4 

1.000 

-5.81 

23.5 

-3.48 

35.4 

1.000 

dawsonite 

— 

— 

— 

-7.00 

62.8 

— 

— 

— 

edolomite 

-3.19 

36.1 

0.500 

-7.53 

52.2 

-5.11 

34.8 

0.500 

fdolomite 

-3.76 

56.7 

0.500 

-8.60 

95.3 

-5.37 

45.7 

0.500 

magnesite 

-6.38 

14.4 

1.000 

-9.34 

23.5 

-5.22 

62.8 

1.000 

a.  Rate  constant  k  computed  from  A  and  E,  25°C,  pH  =  0,  mole  m'2  s'1. 

b.  Arrhenius  activation  energy  E,  kJ  mole'1. 

c.  Reaction  order  n  with  respect  to  H+. 

d.  Reaction  order  n  with  respect  to  P(C02). 

e.  Sedimentary  (disordered)  dolomite. 

f.  Hydrothermal  (ordered)  dolomite. 


3.9  Sulfates 


Dissolution  rate  data  for  several  sulfate  minerals  were  evaluated:  anhydrite  (Barton  and 
Wilde,  1971;  Bildstein  et  al.,  2001;  Dove  and  Czank,  1995;  Jeschke  and  Dreybrodt,  2002b); 
anglesite  (Dove  and  Czank,  1995);  barite  (Christy  and  Putnis,  1993;  Dove  and  Czank,  1995; 
Komicker  et  al.,  1991);  celestite  (Dove  and  Czank,  1995);  and  gypsum  (Jeschke  et  al.,  2001; 
Raines  and  Dewers,  1997).  Data  for  the  sulfate  minerals  were  selected  as  follows:  anhydrite, 
from  Barton  and  Wilde  (1971)  and  Dove  and  Czank  (1995);  anglesite,  barite  and  celestite,  from 
Dove  and  Czank  (1995);  and  gypsum,  from  Raines  and  Dewers  (1997).  The  rate  parameters 
were  regressed  using  eqn.  (11),  and  are  summarized  in  Table  34.  Date  were  insufficient  to 
calculate  activation  energies  for  the  neutral  mechanism  for  gypsum,  and  the  acidic  mechanisms 
for  anglesite  and  barite;  activation  energies  for  the  neutral  mechanisms  for  anglesite  and  barite 
are  assigned  to  the  respective  acidic  mechanisms.  The  acidic  mechanism  for  celestite  apparently 
extends  to  alkaline  pH,  and  there  is  no  evident  neutral  mechanism. 
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Table  34.  Sulfate  Mineral  Dissolution  Rate  Parameters. 


Acid  Mechanism 

Neutral  Mechanism 

alog  k 

bE 

cn 

alog  k 

bE 

anglesite 

-5.58 

31.3 

0.298 

-6.50 

31.3 

anhydrite 

-- 

- 

- 

-3.19 

14.3 

gypsum 

- 

- 

- 

-2.79 

0 

barite 

-6.90 

30.8 

0.220 

-7.90 

30.8 

celestite 

-5.66 

23.8 

0.109 

- 

~ 

a.  Rate  constant  k  computed  from  A  and  E,  25°C,  pH  =  0,  mole  m"2  s'1. 

b.  Arrhenius  activation  energy  E,  kJ  mole'1. 

c.  Reaction  order  n  with  respect  to  H+. 


3.10  Sulfides 


Dissolution  rate  data  for  several  sulfide  minerals  were  evaluated:  AsS,  amorphous 
(Lengke  and  Tempel,  2003);  AS2S3,  amorphous  (Lengke  and  Tempel,  2001);  Orpiment  and 
AS2S3,  amorphous  (Lengke  and  Tempel,  2002);  pyrite  (Gibbs  et  ah,  1997a;  Gibbs  et  al.,  1997b; 
Holmes  and  Crundwell,  2000;  Kamei  and  Ohmoto,  2000;  Lennie  and  Vaughan,  1992;  McKibben 
and  Bames,  1986;  Moses  and  Hennan,  1991;  Moses  et  al.,  1987;  Nicholson  et  al.,  1988; 
Nicholson  et  al.,  1990;  Rickard,  1997;  Rickard  and  Luther,  1997;  Rickard,  1975;  Rimstidt  and 
Vaughan,  2003;  Schoonen  and  Barnes,  1991a;  Schoonen  and  Barnes,  1991b;  Schoonen  and 
Barnes,  1991c;  Wiersma  and  Rimstidt,  1984;  Williamson  and  Rimstidt,  1994);  pyrrhotite 
(Dekkers  and  Schoonen,  1994;  Janzen  et  al.,  2000;  Thomas  et  al.,  2001);  Realgar  (Lengke  and 
Tempel,  2003);  and  Sphalerite  (Weisener  et  al.,  2003). 

Data  for  the  sulfide  minerals  were  selected  as  follows:  AS2S3,  amorphous,  from  Lengke 
and  Tempel  (2001);  pyrite,  from  McKibben  and  Bames  (1986);  and  pyrrhotite,  from  Janzen  et  al. 
(2000).  The  rate  parameters  were  regressed  using  eqn.  (1 1)  or  compiled  directly  from  the 
references,  and  are  summarized  in  Table  35. 


Table  35.  Sulfide  Mineral  Dissolution  Rate  Parameters. 


Acid  Mechanism 

Neutral  Mechanism 

Base  Mechanism 

alog  k 

bE 

cn 

dn 

alog  k 

bE 

en 

alog  k 

bE  cn 

pyrite 

-7.52 

56.9 

-0.500 

0.500 

-4.55 

56.9 

0.500 

- 

- 

fpyrrhotite 

-8.04 

50.8 

-0.597 

0.355 

- 

- 

- 

-- 

- 

9pyrrhotite 

-6.79 

63.0 

-0.090 

0.356 

- 

- 

- 

-- 

-- 

hAs2S3 

- 

- 

- 

- 

-9.83 

8.7 

0.180 

-17.39 

8.7  -1.208 

a.  Rate  constant  k  computed  from  A  and  E,  25°C,  pH  =  0,  mole  m'2  s'1. 

b.  Arrhenius  activation  energy  E,  kJ  mole'1. 

c.  Reaction  order  n  with  respect  to  H+. 

d.  Reaction  order  n  with  respect  to  Fe+++. 

e.  Reaction  order  n  with  respect  to  02. 

f.  Monoclinic  pyrrhotite. 

g.  Hexagonal  pyrrhotite. 

h.  Amorphous  As2S3. 
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3.11  Phosphates 


Dissolution  rate  data  for  apatite  from  several  references  were  evaluated  (Chien,  1977; 
Chien  et  al.,  1980;  Christoffersen  et  al.,  1978;  Fawzi  et  al.,  1978;  Fox  et  al.,  1978;  Guidry  and 
Mackenzie,  2003;  Hull  and  Hull,  1987;  Valsami- Jones  et  al,  1998).  Data  for  hydroxyapatite  and 
fluorapaite  were  selected  from  Fawzi  et  al.  (1978)  and  Valsami-Jones  et  al.  (1998).  The  rate 
parameters  were  regressed  using  eqn.  (11),  and  are  summarized  in  Table  36.  For  the  neutral 
mechanism,  a  “cut-off’  rate  constant  at  25  °C  was  calculated  so  that  rates  do  not  decrease  with 
pH  increasing  above  7,  and  the  activation  energy  is  set  equal  to  the  acidic  mechanism.  The 
uncertainty  in  activation  energy  may  be  quite  large  because  it  was  computed  from  dissolution 
rates  obtained  from  experiments  with  minimal  temperature  difference,  25  and  28  °C  (Fawzi  et 
al.,  1978);  these  data  were  for  experiments  with  hydroxyapatite,  and  the  tabulated  activation 
energy  for  fluorapatite  dissolution  is  taken  equal  to  that  of  hydroxyapatite. 


Table  36.  Phosphate  Mineral  Dissolution  Rate  Parameters. 


Acid  Mechanism 

Neutral  Mechanism 

alog  k 

bE 

cn 

alog  k 

bE 

hydroxyapatite 

-4.29 

250.0 

0.171 

-6.00 

250.0 

fluorapatite 

-3.73 

250.0 

0.613 

-8.00 

250.0 

a.  Rate  constant  k  computed  from  A  and  E,  25°C,  pH  =  0,  mole  m'2  s'1. 

b.  Arrhenius  activation  energy  E,  kJ  mole"1. 

c.  Reaction  order  n  with  respect  to  H+. 
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3.12  Halides 


Dissolution  rate  data  for  halite  (Alkattan  et  al.,  1997a;  Alkattan  et  al.,  1997b)  and  fluorite 
(Bosbach  et  al.,  1995;  Christoffersen  et  al.,  1988;  Gardner  and  Nancollas,  1976;  Hamza  et  al., 
1987;  Hamza  and  Hamdona,  1991;  Zhang  and  Nancollas,  1990)  were  evaluated.  Data  for  halite 
were  selected  from  Alkattan  (1997a),  and  the  rate  parameters  were  regressed  using  eqn.  (11);  the 
far  from  equilibrium  rate  constant  was  calculated  from  data  at  f2  =  0.90  to  0.94  and  the  apparent 
rate  constant.  For  fluorite,  the  rate  constant  at  25  °C  in  pure  water  was  selected  from  Zhang  and 
Nancollas  (1990),  where  the  far  from  equilibrium  rate  constant  was  calculated  from  data  at  D  = 
0.40  and  the  apparent  rate  constant,  and  the  activation  energy  was  selected  from  Hamza  and 
Hamdone  (1991).  The  rate  parameters  for  halite  and  fluorite  are  summarized  in  Table  37. 

Table  37.  Halide  Mineral  Dissolution  Rate  Parameters. 


Acid  Mechanism 

Neutral  Mechanism 

alog  k 

bE  cn 

alog  k 

bE 

halite 

- 

.. 

-0.21 

7.4 

fluorite 

-7.14 

73.0  1.000 

-13.79 

73.0 

a.  Rate  constant  k  computed  from  A  and  E,  25°C,  pH  =  0,  mole  rrf2  s'1. 

b.  Arrhenius  activation  energy  E,  kJ  mole  f 

c.  Reaction  order  n  with  respect  to  H+. 
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4.  CONCLUSIONS 


This  work  is  intended  as  a  preliminary  starting  point  for  the  implementation  of  the 
kinetics  of  mineral  dissolution  and  precipitation  in  aqueous  fluids  in  reaction  path  modeling  of 
water-rock-gas  systems.  The  list  of  cited  references  may  be  as  useful  as,  or  perhaps  more  useful 
than,  the  rate  parameters  compiled  herein,  as  the  references  commonly  provide  infonnation  to 
further  refine  rate  equations  that  have  been  implemented,  particularly  with  regard  to  the  use  of 
adsorption  isotherms  to  quantify  the  effects  of  dissolved  constituents  upon  reaction  rates.  The 
general  rate  equation  and  rate  parameters  were  compiled  specifically  for  the  program 
GAMSPATH  (Perkins  et  ah,  1997),  which  simulates  an  infinitely  well  stirred  batch  reactor.  For 
the  program  GAMSPATH,  the  rate  parameters  compiled  provide  the  ability  to  compute,  to  a  first 
approximation,  the  minimum  length  of  time  required  for  a  system  to  equilibrate. 

For  most  minerals,  there  is  currently  a  high  degree  of  uncertainty  in  experimentally 
determined  dissolution  and  precipitation  rates.  For  many  minerals,  variation  in  composition, 
degree  of  cation  disorder,  degree  of  crystallinity,  as  well  as  the  frequency  and  distribution  of 
crystal  defects,  contribute  to  these  uncertainties.  In  many  cases,  the  lack  of  data  over  a  wide 
range  of  temperatures  results  in  a  large  degree  of  uncertainty  in  the  Arrhenius  activation  energy. 
Furthermore,  the  possible  effects  of  many  common  dissolved  species  have  not  been  quantified. 

It  is  clear  that  much  experimental  work  remains  to  be  done  to  detennine  the  kinetics  of  mineral 
dissolution  and  precipitation  under  conditions  that  are  encountered  in  natural  systems. 
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